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THE COMPRESSIBLE VISCOUS SURFACE-INTERNAL WAVE PROBLEM: 
STABILITY AND VANISHING SURFACE TENSION LIMIT 

JUHI JANG, IAN TICE, AND YANJIN WANG 


Abstract. This paper concerns the dynamics of two layers of compressible, barotropic, viscous 
fluid lying atop one another. The lower fluid is bounded below by a rigid bottom, and the upper 
fluid is bounded above by a trivial fluid of constant pressure. This is a free boundary problem: 
the interfaces between the fluids and above the upper fluid are free to move. The fluids are 
acted on by gravity in the bulk, and at the free interfaces we consider both the case of surface 
tension and the case of no surface forces. We establish a sharp nonlinear global-in-time stability 
criterion and give the explicit decay rates to the equilibrium. When the upper fluid is heavier 
than the lower fluid along the equilibrium interface, we characterize the set of surface tension 
values in which the equilibrium is nonlinearly stable. Remarkably, this set is non-empty, i.e. 
sufficiently large surface tension can prevent the onset of the Rayleigh-Taylor instability. When 
the lower fluid is heavier than the upper fluid, we show that the equilibrium is stable for all 
non-negative surface tensions and we establish the zero surface tension limit. 


1. Introduction 

1.1. Formulation in Eulerian coordinates. We consider two distinct, immiscible, viscous, 
compressible, barotropic fluids evolving in a moving domain Q(t) = _(t) U 0_(i) for time 

t > 0. One fluid (+), called the “upper fluid”, fills the upper domain 

M+(t) = {y € T 2 x M | rj-(yi,y 2 ,t ) < y 3 < l + rj + (yi, y 2 , t)}, (1.1) 

and the other fluid (—), called the “lower fluid”, fills the lower domain 

= {y € T 2 x R | -b < y 3 < r)-(yi, y 2 , t)}. (1.2) 

Here we assume the domains are horizontally periodic by setting T 2 = ( 2ttLiT ) x ( 2nL 2 T) for 
T = M/Z the usual 1-torus and L±,L 2 > 0 the periodicity lengths. We assume that £, b > 0 
are two fixed and given constants, but the two surface functions y± are free and unknown. The 
surface T + (f) = {y 3 = l + r] + (yi,y 2 ,t)} is the moving upper boundary of H_|_(i) where the upper 
fluid is in contact with the atmosphere, T_(t) = {y 3 = rj_(yi,y 2l t)} is the moving internal 
interface between the two fluids, and Ej, = {y 3 = —b} is the fixed lower boundary of fi_(t). 

The two fluids are described by their density and velocity functions, which are given for each 
t > 0 by p±(-,t) : fl-i-(f) -A M + and u±(-,t) : Q±(t) -A M 3 , respectively. In each fluid the 
pressure is a function of density: P± = P±(p±) > 0, and the pressure function is assumed to be 
smooth, positive, and strictly increasing. For a vector function u € M 3 we define the symmetric 
gradient by ( Ou)ij = diUj + djUi for i,j = 1,2,3; its deviatoric (trace-free) part is then 

2 

D°u = Du — - div ul, (1.3) 

3 

where I is the 3x3 identity matrix. The viscous stress tensor in each fluid is then given by 

§±(h±) := / u-j-lD) 0 ft-|- + p,'± div u±I, (1-4) 

where p± is the shear viscosity and p!± is the bulk viscosity; we assume these satisfy the usual 
physical conditions 

P± > 0, p!± > 0. (1.5) 
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The tensor P±(p±)I — §±(^±) is known as the stress tensor. The divergence of a symmetric 
tensor M is defined to be the vector with components (divM)j = Note then that 

div (P±(p±)I — S±(fi±)) = VP±(p ± ) — p±Au± - + //±) V div u±. (1.6) 

For each t > 0 we require that (p±, u±, p±) satisfy the following equations 


' d t p ± + div(p±u±) = 0 

p±{d t u± + u± • Vu±) + VP±(p ± ) - div§±(tt±) = -gp±e 3 
dtP± = u 3) ± - ui,±d yi r]± - u 2 ,±dy 2 r]± 

< ( P+(p+)I - § + (u+))n + = pat m n + - a + n + n + 

(P + (p + )I — §+(?2 + ))n_ = (P-(p-)I — §_(n_))n_ + 
u .|_ = u- 
,u- = 0 


in Q±(t) 
in 

on T±(t) 
on T+(t) 
on T_(t) 
on T_(t) 
on S;,. 


(1.7) 


In the equations —gp±e 3 is the gravitational force with the constant g > 0 the acceleration of 
gravity and e 3 the vertical unit vector. The constant p a tm > 0 is the atmospheric pressure, 
and we take cr± > 0 to be the constant coefficients of surface tension. In this paper, we let 
V* denote the horizontal gradient, div* denote the horizontal divergence and A* denote the 
horizontal Laplace operator. Then the upward-pointing unit normal of T-j-(t), n±, is given by 


n± 


( —V*7?±,l) 
V 7 ! + |V*?7±| 2 


( 1 . 8 ) 


and T~L±, twice the mean curvature of the surface T±(t), is given by the formula 


U± 


div* 


( V*?7± \ 

V\/l + |V*r 7 ±| 2 J ' 


(1.9) 


The third equation in (11.711 is called the kinematic boundary condition since it implies that the 
free surfaces are advected with the fluids. The boundary equations in (11.71) involving the stress 
tensor are called the dynamic boundary conditions. Notice that on T_(f), the continuity of 
velocity, u+ = v,-, means that it is the common value of u± that advects the interface. For 
a more physical description of the equations and the boundary conditions in (11.71) . we refer to 


To complete the statement of the problem, we must specify the initial conditions. We suppose 
that the initial surfaces r -i- (0) are given by the graphs of the functions ?7±(0), which yield the 
open sets H±(0) on which we specify the initial data for the density, /5±(0) : fl±(0) —>• M + , and 
the velocity, ft±(0) : fl±(0) -A- M 3 . We will assume that t + p + ( 0) > r]-( 0) > —b on T 2 , which 
means that at the initial time the boundaries do not intersect with each other. 


1.2. Equilibria. We seek a steady-state equilibrium solution to CL 2 D with u± = 0,r/± 
the equilibrium domains given by 

12+ = {y € T 2 x R | 0 < 2/3 < 1} and = {y € T 2 x R | — b < y 3 < 0}. 


Then (11.71) reduces to an ODE for the equilibrium densities p± = p±(y 3 )'- 


d(P+(p+)) 

dy 3 

d(P~(p-)) 


= -9 P+, 


= -9P-, 


dy 3 

P+iP+i^Y) = Patmi 

[p+ip+m = p-(p-(o)). 


for y 3 € ( 0 , £), 
for y 3 € (- 6 , 0 ), 


0 , and 

( 1 . 10 ) 


( 1 . 11 ) 


The system (11.111) admits a solution p± > 0 if and only if we assume that the equilibrium 
heights b,£ > 0, the pressure laws P±, and the atmospheric pressure p a tm satisfy a collection of 
admissibility conditions, which we enumerate below. 
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First, in order for the third equation in (11.111) to be achievable, we must have that p a tm € 
P+(R + ). This then allows us to uniquely define 

Pi := P+ l (Patm) > 0. (1.12) 

We then introduce the enthalpy function h .|_ : (0, oo) —» M given by 


h+(z) = 


f- 

J p\ 


P' + {r) 


dr , 


(1.13) 


which is smooth, strictly increasing, and positive on (pi,oo). The first ODE in (11.111) then 
reduces to dh + (p + )/dy 3 = —g , subject to the condition that p+(£) = p\ , and hence 

P+{y 3) = h^-(h + (pi) + g(l - 2/3)) = h+ l (g{£ - y 3 )). (1.14) 

This gives a well-defined, smooth, and decreasing function p + : [0,6] —>• [, 01 , 00 ) if and only if 


1 poo pr ( r \ 

gt £ /i + ((pi,oo)) <(=)> 0 < t < - / ——-dr. 

W* r 


(1.15) 


The equation (11.141) yields the value of the upper density at the equilibrium internal interface: 

p + := p+(0) = h~ l {gl). (1.16) 

Then the last equation in ([1.11D is achievable if and only if P + (p + ( 0)) = P + (h'^ 1 (g£)) € P_( 


and in this case yields the value of the lower density at the equilibrium internal interface: 

p~ := p_(0) = Pl\P + (p + )) = Pl\h- + \gl)). (1.17) 

Then we define /i_ : (0, 00 ) —> M via 


h _( Z )=r^±dr 


(1.18) 

(1.19) 


and solve the ODE for p-(y 3 ) as before to obtain 

P-{y 3 ) = hz 1 (h-(p~) - gy 3 ) = hZ l {-gy 3 ). 

This gives a well-defined, smooth, and decreasing function : [—6,0] —>• [p _ ,oo) if and only if 

f *°° P' (r) 


1 r 

gb € h-((p , 00 )) <(=)> 0 < 6 < - 

9 Js- 


-dr. 


( 1 . 20 ) 


0 <£< 


-dr, 


( 1 . 21 ) 


For the sake of clarity we now summarize the admissibility conditions that are necessary and 
sufficient for the existence of an equilibrium: 

(1) Patm € -P+(M + ), which defines pi := P^iPatm)- 

(2) t satisfies 

1 P[{r) 

9 Jpi 

which defines p + := h^ 1 (g£). 

(3) £ also satisfies P + (h+ l (g£)) € i 3 _(M + ), which defines p~ := PI 1 (h^ 1 (g£)). 

(4) 6 satisfies 

1 PZ{r) , 

0 < 6 < - / - -dr. 

9 Jp- r 

As an example, consider the case 

P±(z) = K±z a± for K± > 0 and a± > 1. 

Then for any p a tm ,6, 6 > 0 we may solve to see that 

fpatm\ {a+ ~ 1)/a+ , g(a+-l)y_A 


( 1 . 22 ) 

(1.23) 


P+(ya) = 

P-{ys) = 


\K+J 


+ 


K + a + 


K + (p + m°+\ fa - iy °- , g(a-~l) ( 


K_ 


1 !/(«--!) 


AT_a_ 


(1.24) 
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are the desired solutions to (11.111) . Note in particular that in this case the admissibility condi¬ 
tions are satisfied trivially. 

Throughout the rest of the paper we will assume that the above admissibility conditions are 
satisfied for a fixed £, b > 0, which then uniquely determine the equilibrium density p±(y 3 ). In 
turn these fix the equilibrium masses 

r£ Att 2 LiL2 


M, = 


M_ = 


/ P+ = (4tt 2 LiL 2 ) / p + (y 3 )dy 3 - 
Jq + Jo 

f p_ = (4vr 2 LiL 2 ) [ p-{y 3 )dy 3 
J-b 


{P+{h + l (g£)) -Patm), 


9 

47t 2 LiL 2 
9 


(1.25) 


{P-{hZ\gb)) - P + {h+\gl))) 


It is perhaps more natural to specify M + ,M_ > 0 and then to determine £, b > 0 in terms of 
M + , M_,p a tm and the pressure laws by way of (11.251) . This can be done under the assumption 
that M_|_, M_, Patm satisfy certain admissibility conditions similar to those written above for 
£, 5, p a tm- For the sake of brevity we will neglect to write these explicitly. 

Notice that the equilibrium density can jump across the internal interface. The jump, which 
we denote by 

M '■= P+{0) ~ P-(Q) = P + ~ P~, (1-26) 

is of fundamental importance in the the analysis of solutions to (ED near equilibrium. Indeed, 
if M > 0 then the upper fluid is heavier than the lower fluid along the equilibrium interface, 
and the fluid is susceptible to the well-known Rayleigh-Taylor gravitational instability. 


1.3. History and known results. Free boundary problems in fluid mechanics have attracted 
much interest in the mathematical community. A thorough survey of the literature would prove 
impossible here, so we will primarily mention the work most relevant to our present setting, 
namely that related to layers of viscous fluid. We refer to the review of Shibata and Shimizu 
[ 22 ] for a more proper survey of the literature. 

The dynamics of a single layer of viscous incompressible fluid lying above a rigid bottom, 
i.e. the incompressible viscous surface wave problem, have attracted the attention of many 
mathematicians since the pioneering work of Beale [3]. For the case without surface tension, 
Beale [3] proved the local well-posedness in the Sobolev spaces. Hataya [12] obtained the global 
existence of small, horizontally periodic solutions with an algebraic decay rate in time. Guo 
and Tice El El dD] developed a two-tier energy method to prove global well-posedness and 
decay of this problem. They proved that if the free boundary is horizontally infinite, then the 
solution decays to equilibrium at an algebraic rate; on the other hand, if the free boundary is 
horizontally periodic, then the solution decays at an almost exponential rate. The proofs were 
subsequently refined by the work of Wu [35]. For the case with surface tension, Beale [3j proved 
global well-posedness of the problem, while Allain [1] obtained a local existence theorem in two 
dimensions using a different method. Bae [2] showed the global solvability in Sobolev spaces via 
energy methods. Beale and Nishida [5] showed that the solution obtained in [3] decays in time 
with an optimal algebraic decay rate. Nishida, Teramoto and Yoshihara [19] showed the global 
existence of periodic solutions with an exponential decay rate in the case of a domain with a 
flat fixed lower boundary. Tani [28] and Tani and Tanaka [29] also discussed the solvability 
of the problem with or without surface tension by using methods developed by Solonnikov in 
[23 . 24; |25]. Tan and Wang [26] studied the vanishing surface tension limit of the problem. 

There are fewer results on two-phase incompressible problems, i.e. the incompressible viscous 
surface-internal wave or internal wave problems. Hataya m proved an existence result for a 
periodic free interface problem with surface tension, perturbed around Couette flow; he showed 
the local existence of small solution for any physical constants, and the existence of exponentially 
decaying small solution if the viscosities of the two fluids are sufficiently large and their difference 
is small. Priiss and Simonett [20] proved the local well-posedness of a free interface problem 
with surface tension in which two layers of viscous fluids fill the whole space and are separated 
by a horizontal interface. For two horizontal fluids of finite depth with surface tension, Xu and 
Zhang [36] proved the local solvability for general data and global solvability for data near the 
equilibrium state using Tani and Tanaka’s method. Wang and Tice [32] and Wang, Tice and 
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Kim [33] adapted the two-tier energy methods of la Gama to develop the nonlinear Rayleigh- 
Taylor instability theory for the problem, proving the existence of a sharp stability criterion 
given in terms of the surface tension coefficient, gravity, periodicity lengths, and [pj. 

The free boundary problems corresponding to a single horizontally periodic layer of com¬ 
pressible viscous fluid with surface tension have been studied by several authors. Jin m and 
Jin-Padula [16] produced global-in-time solutions using Lagrangian coordinates, and Tanaka 
and Tani m produced global solutions with temperature dependence. However, to the best of 
our knowledge, the problem for two layers of compressible viscous fluids remains open except 
for the linear analysis of Guo and Tice [7]. 

The two-layer problem is important because it allows for the development of the classical 
Rayleigh-Taylor instability EH EDI, at least when the equilibrium has a heavier fluid on top 
and a lighter one below and there is a downward gravitational force. In this paper we identify 
a stability criterion and show that global solutions exist and decay to equilibrium when the 
criterion is met. In our companion paper m we show that the stability criterion is sharp, 
as in the incompressible case [321 [33] . and that the Rayleigh-Taylor instability persists at the 
nonlinear level (the linear analysis was developed in [7]). 

2. Reformulation and main results 

In this section, we reformulate the system m by using a special flattening coordinate 
transformation and deliver the precise statements of the main results. 

2.1. Reformulation in flattened coordinates. The movement of the free surfaces T±(t) 
and the subsequent change of the domains Q±(t) create numerous mathematical difficulties. To 
circumvent these, we will switch to coordinates in which the boundaries and the domains stay 
fixed in time. Since we are interested in the nonlinear stability of the equilibrium state, we 
will use the equilibrium domain. We will not use a Lagrangian coordinate transformation, but 
rather utilize a special flattening coordinate transformation motivated by Beale [1]. 

To this end, we define the fixed domain 

H = fl + U with fl_|_ := {0 < X 3 < 1 } and := {—b < X 3 < 0 }, ( 2 - 1 ) 

for which we have written the coordinates as x € Cl. We shall write := {£3 = £} for the 
upper boundary, £_ := {£3 = 0} for the internal interface and := {x^ = —6} for the lower 
boundary. Throughout the paper we will write £ = £_|_ U £_. We think of r]± as a function on 
£± according to r] + : (T 2 x {£}) x R + —> R and r: (T 2 x {0}) x R + -A R, respectively. We 
will transform the free boundary problem in Cl(t) to one in the fixed domain Cl by using the 
unknown free surface functions T]±. For this we define 

f) + := V+rj + = Poisson extension of r] + into T 2 x {^3 < £} (2-2) 

and 

fj_ \= V-T= specialized Poisson extension of r/_ into T 2 x R, (2-3) 

where V± are defined by (IA.4I) and HA.91) . Motivated by Beale [3], the Poisson extensions fj± 
allow us to flatten the coordinate domains via the following special coordinate transformation: 

(x 1 ,x 2 ,x 3 + bifj + + b 2 fj-) := @(t,x) = ( 2 / 1 , 2 / 2 , 2 / 3 ) € fi±(i), (2.4) 

where we have chosen b\ = bi{x^,) 1 b 2 = b 2 {x^) to be two smooth functions in R that satisfy 

6i(0) = &!(-&) = 0, bi(£) = 1 and b 2 {l) = b 2 {-b) = 0,6 2 (0) = 1. (2.5) 

Note that ©(£+,t) = r+(i), ©(£_,£) = r_(t) and @(-,2) |e 6 = Id |s 6 . 

Note that if rj is sufficiently small (in an appropriate Sobolev space), by Lemma fA.il and usual 
Sobolev embeddings, then the mapping 0 is a diffeomorphism. This allows us to transform the 
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problem (11.71) to one in the fixed spatial domain 17 for each t > 0. In order to write down the 
equations in the new coordinate system, we compute 

/ 1 0 0 \ / i ° -AK \ 

V0 = I 0 1 0 and A := (V0 _1 ) T = I 0 1 -BK . (2.6) 

\ A B J J \ 0 0 K ) 

Here the components in the matrix are 

A = d^e, B = d 2 9, J = 1 + d 3 0 , K = J -1 , (2.7) 

where we have written 

9 := bifj+ + b 2 fj— ( 2 . 8 ) 

Notice that J = detV0 is the Jacobian of the coordinate transformation. It is straightforward 
to check that, because of how we have defined p- and 0, the matrix A is regular across the 
interface £_. 

We now define the density p± and the velocity u± on by the compositions p±(x,t ) = 
p±(@±(x,t),t) and u±(x,t) = u±(Q±(x,t),t)- Since the domains 17-t and the boundaries 
are now fixed, we henceforth consolidate notation by writing / to refer to f± except when 
necessary to distinguish the two; when we write an equation for / we assume that the equation 
holds with the subscripts added on the domains or To write the jump conditions on 
E_, for a quantity f = f±, we define the interfacial jump as 

lfj:=f + \ { x 3 = 0 }-f-\ {x3=0} . (2.9) 

Then in the new coordinates, the PDE (11.71) becomes the following system for ( p,u,p ): 


'dtp — Kdt9d 3 p + div^ (pu) = 0 in 17 

p(d t u - Kd t 9d 3 u + u ■ V^u) + V A P(p) - div^§ A (u) = -gpe 3 in 17 
dtp = u • Af on E 

< ( P(p)I - Sa(u)W = PatmAf - cr+TLAT on S + 

\P(p)I — §a( w )1 Af = a-TLAT on S_ 

[it] = 0 on S_ 

^U- = 0 on Sb. 


( 2 . 10 ) 


Here we have written the differential operators V_ 4 , div. 4 , and A _4 with their actions given by 
(V A f)i ■■= Aijdjf, div _4 X := AijdjXt, and A A f := div^ V A f ( 2 . 11 ) 

for appropriate / and X. We have also written 

Af := {~d\p, ~d 2 p, 1 ) ( 2 . 12 ) 

for the non-unit normal to E(t), and we have written 


S A ,±( u ) := P-±^ ) A U + p'± divA ul, = B>au — — div_4 ul , 

O 

and (ID^u)^' — Atkd^Uj T Ajkdj^ut. (2.13) 

Note that if we extend div^ to act on symmetric tensors in the natural way, then div .4 = 
pA A u + (p/3 + p')V a div ^4 u. Recall that A is determined by p through (12.61) . This means that 
all of the differential operators in ( 12 . 101 ) are connected to p , and hence to the geometry of the 
free surfaces. 


Remark 2.1. The equilibrium state given by (11.111) corresponds to the static solution ( p,u,p ) 

(p,o,o) of dZrnD- 
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2 . 2 . Local well-posedness. As mentioned in the introduction, the purpose of this article 
is to establish a sharp nonlinear stability criterion for the equilibrium state (p, 0 , 0 ) in the 
compressible viscous surface-internal wave problem ( 12 . 101 ) . Since the equilibrium velocity and 
free surfaces are trivial, it suffices to study u and r/ directly. However, since the equilibrium 
density is nontrivial, we must define a perturbation of it. We do so by defining 

q:=p-p-d 3 p9. (2.14) 

Note that if we consider the natural density perturbation g = p — p, then the resulting PDE 
has a linear term —gge 3 in the right hand side of the second momentum equation, which is not 
convenient for employing elliptic regularity theory; our special density perturbation q defined 
by (j2.14l) eliminates this linear term (see (|3.3p ). For the sake of brevity, we will not record here 
the equations for (g, u, rj) obtained from plugging this perturbation into ( 12 . 101 ) : they may be 
found later in the paper in ()3.24l) . 

Before discussing our main results we mention the local well-posedness theory for ()2.101) . 
perturbed around the equilibrium state. The energy space in which local solutions exist is 
defined in terms of certain functionals, which are sums of Sobolev norms. We refer to Section 
1231 for details on our notation for Sobolev spaces and norms. For a generic integer n > 3, we 
define the energy as 


£n := E 

3=0 


+ 


lL 


+ E II dt q 

3 = 1 


2n—2j 

+ 'H([pJ)mm{l,cr + ,cr- - cr c } 
and the corresponding dissipation as 


2n-2j+l 


Elh 

3= 1 


2n-2j+3/2 


l2n+l 


+ %(— [/?]) 


12 n 


+ min{l, a} 


12n+l 


(2.15) 


V a ■= 


E 1%* 

3= 0 


+ 


71+1 

2n—2j+l ML-i+ElK? 

3 =2 
2 


lL 


2n—2j+2 

2i II„ll2 


+ W(H)min{l,<T + ,<T_ -a c ,a + , (<r_ - a c ) } |M | 2n+3/2 
+ ft(- IpD (lMlL-1/2 + min + O’ 2 } lkllL+3/2) 


n +1 


+ \\dtV\\ln-l/2 + V 2 \\dtV\\ln+l/2 + E 

3 =2 


2n—2j+5/2 


, (2.16) 


where here T~L = X(o,oo) is the Heaviside function. Throughout the paper we will consider both 
n = 2 N and n = N + 2 for the integer N > 3. We also define 

?2N-.= h\\l N+ l/2- (2-17) 

Note that we have included the surface tension coefficient a± in the definitions (|2.151) and (12.161) 
so that we will be able to treat the cases with and without surface tension together. Throughout 
the paper, we assume that cr± have the same sign: either cr + = cr_ = 0 or else a± > 0. Note 
that this assumption excludes the cases of cr_ = 0, > 0 and <r_ > 0, cr + = 0. 

In order to produce local solutions on [0,T] for which sup 0<t<T £% N (t) remains bounded, we 
must assume that the initial data (go, uo, t?o) satisfy 2N systems of compatibility conditions (a 
system for each time derivative of order 0 to 2N — 1). These are natural for solutions to (12.101) 
in our energy space, but they are cumbersome to write. We shall neglect to record them in this 
paper and instead refer to our companion paper m for their precise enumeration. 

Our local existence results for (12.101) is then as follows. 


Theorem 2.2. Assume that either a± = 0 or else a± > 0. Let N >3 be an integer. Assume 
that the initial data (g 0 ,«o,i?o) satisfy HuollLv + IMlLv + II+^at+i^ + G II ^*Vo\\In < 00 as 
well as the compatibility conditions enumerated in |13j. There exist 0 < <5o, To < 1 so that if 
H+lLv + 11 go 11 47 v + ll^/o II 47 V+ 1/2 + a HV^ollLv < 5o and 0 < T < T 0; then there exists a unique 
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triple ( q,u,r /), achieving the initial data, so that (p = p + q + d 3 pO, u, p) solve (12.101) . The 
solution obeys the estimates 


sup £% N (t) + sup J~ 2 n(£) + 
0 <t<T 0 <t<T 



p(t)J(t)dr+ l u(t) 


2 

dt 

(0 H 1 )* 


< 


K(cr, [p]) (lluoll 47V + 1 1 Q'0 1147V + II^O ||liV+l/2 + a 


i2 

I4JV 


where K (<7, [pj) is a constant of the form (12.401) . 


(2.18) 


Theorem 12.21 can be deduced readily from Theorem 2.1 of m- Indeed, Theorem 2.1 of m 
is stated in more general form, where we only require ||?7o II4JV—1/2 to be small, and no smallness 
condition is imposed on uq or go- We record the version of local well-posedness in Theorem 12.21 
so that it can be employed directly in proving our global well-posedness result. 


2.3. Main results. To state our results, we now define some quantities. For a given jump 
value in the equilibrium density [[/?], we first define the critical surface tension value by 

ch := [p] gmax{Lf,L|}, (2.19) 

where we recall that L\, L 2 are the periodicity lengths. Note that this is the same critical value 
identified for the incompressible problem in [32] [33 ]. 

For the global well-posedenss theory, we assume further that the initial masses are equal to 
the equilibrium masses: 

/ Po,+Jq,+ = / P+ and / po,_J 0 ,- = / p— (2.20) 

J f2_j_ J Q + J J Q_ 

Note that for sufficiently regular solutions to the problem (12.101) (and in particular our global 
solutions), the condition (12.201) persists in time, i.e. 


[ p+(t)J+(t) = I p+ and [ p-{t)J-{t) = [ P- 
J r2_i_ J j £i— 


( 2 . 21 ) 


Indeed, from the first continuity equation and the equations dtp = u ■ J\f on X and u_ = 0 on 
Xf,, we find, after integrating by parts, that 

d f P+J+= f p+d t J++ d t p+J+= j p + d t d 3 0 + + d t d + d 3 p + + J + div A+ (p + u + ) 

J 0+ 4 o_(_ ^ 22^ 

= / P+(dfn+ - U + ■ A4) - / P+{dtV- ~ U ■ M-) = 0, 

Jt,+ Jt,- 


dt 


and 




P- (dtP- - u-N -) = 0. 


(2.23) 


The conservation of masses (12.211) plays a key role when [p]] >0, as it allows us to show the 
positivity of the energy in the stability regime. Moreover, we are interested in showing p(t) —> 0 
and q(t) —>• 0 (and hence p(t) —> p) as t —> 00 in a strong sense; due to the conservation of 
(12.211) . we cannot expect this unless (12.201) is satisfied at initial time. 

We will employ a functional that combines control of the energy functionals given by ([2.151) . 
the dissipation functionals given by ()2. 16[) , and the auxiliary energy given by (12.171) . We define 

Gwit) : = SU P ^2iv( r )+ / T> 2 N (r)dr+ sup (1 + r) 4JV-8 £ft +2 (r) + sup ^ | . (2.24) 

We can now state our global well-posedness result for (12.101) with or without surface tension. 


Theorem 2.3. Suppose that one of the following three cases holds: 


[pj < 0, cr + = 0, = 0 

(2.25) 

IpI < 0, cr + > 0, > 0 

(2.26) 

IpJ > 0, (T + > 0, <7_ > a c , 

(2.27) 
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where in the last case a c > 0 is defined by (|2.19l) . Let N > 3 be an integer. Assume that the 
initial data ( qo,uo,Vo ) satisfy Sf N (0) + < oo as well as the compatibility conditions of 

Theorem \2.2\ and the mass condition (12.2011 . Let Jo > 0 be the constant from Theorem \2.2\ . 

There exists a constant 0 < k < Jo such that 1 /k = K(a,\p\), where K(a , [p] ) is a constant 
of the form (I2.40|) . such that if Sf N ( 0) + J-2jv(0) < k, then there exists a unique triple ( q,u,r]), 
achieving the initial data, such that (p = p + q + d$pQ,u,rj) solve (12.1011 on the time interval 
[0,oo). The solution obeys the estimate 

Q2 N { oo) < K(a, [p]) (Sf N (0) + J- 2iV (0)). (2.28) 

With surface tension, i.e. cases (12.2611 and (12.271) . there exists a constants K(a,\p\), of the 
form (I2.40p . and M(ct, [p]) 7 defined by (|8.46l) and (J8.47|) . so that 

sup 
t> o 

Theorem 12.7)1 is proved later in the paper in Section [9j Before stating our second main result, 
some remarks are in order. 


exp 


K{cr, \p})M(o, 


S' 


N +2 


(*) 


< 

J 


K(a, [pl)^ +2 (0) 


(2.29) 


Remark 2.4. Without surface tension (case (12.251 ) ). the estimate (|2.28l) guarantees that 


sup(l +t) 4N - 
t> o 


II2A+4 + II Q'(^) II22V+4 


II 7 7(^) II27V+4 


<&(0)+^(0)) <«. 


(2.30) 


Since N may be taken to be arbitrarily large, this decay result can be regarded as an “almost 
exponential” decay rate. 

With surface tension (cases (12.261) and (12.27 [) ), the estimate (I2.29p implies that 


sup 

t>o 


exp 


(i^a, \p])M(a, [p]) 


II2V+4 + II9WII2JV+4 + c ( a ’ M) \\v(t)\\ 2A+5) I 

< K{a, lPm% N (0) + -T 27 v(0)) < K(a, [p])«, 


(2.31) 


where C(o,\p\) = %([p]) min{l, a + , cr_ — cr c } + H{— [p|) min{l, cr+, cr_} andTL = X(0,oo) 
usual Heaviside function. This is exponential decay. We then see that one of the main effects 
of surface tension is that it induces a faster decay to equilibrium. 


Remark 2.5. In our companion paper m we show that if the stability criterion is not satis¬ 
fied, then the local solutions are nonlinearly unstable. This establishes sharp nonlinear stability 
criteria for the equilibrium state in the compressible viscous surface-internal wave problem. We 
summarize these and the rates of decay to equilibrium in the following table. 



M <0 

M =0 

M >0 

a± = 0 

nonlinearly stable 
almost exponential decay 

locally well-posed 

nonlinearly unstable 

0 < <7_|_ 

0 < <r < a c 

nonlinearly stable 
exponential decay 

nonlinearly stable 
exponential decay 

nonlinearly unstable 

0 < <7+ 

a c = am 

nonlinearly stable 
exponential decay 

nonlinearly stable 
exponential decay 

locally well-posed 

0 < a + 
a c < am 

nonlinearly stable 
exponential decay 

nonlinearly stable 
exponential decay 

nonlinearly stable 
exponential decay 


Remark 2.6. If we introduce the critical density jump value by 


Me : = 


<7_ 

g max{L^, L(f\ 


> 0 


(2.32) 


for a given lower surface tension value <r_, then Theorem \2.2\ characterizes the stable equilibria 
in terms of M c ; equilibrium is nonlinear stable if [p] < [p] c . In particular, this immediately 
implies that if [pj < 0, the equilibrium is stable for all surface tension values a± > 0. 
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Remark 2.7. Our methods could be readily applied to the internal-wave problem, i.e. the 
problem (12.1011 posed with a rigid top in place of the upper free surface. In this case Theorem 
1 2. -71 again holds, and we identify a nonlinear stability criterion for the compressible viscous 
internal-wave problem. This agrees with the linear stability criterion identified in [7]. 

Remark 2.8. The high regularity context of Theorem \2.3\ is not entirely necessary in the case 
with surface tension. Indeed, the theorem shows in this case that the integer N plays no signif¬ 
icant role in determining the decay rate, which is exponential. We believe that we could refine 
the analysis in the case with surface tension to match the lower regularity context used in the in¬ 
compressible case |32j |33]. We have chosen to forgo this approach here for two reasons. First, it 
allows us to analyze the case with surface tension and the case without simultaneously. Second, 
it allows us to establish the vanishing surface tension limit. 

In our analysis leading to the proof of Theorem 12.31 we have made a serious effort to track 
the dependence of various constants on the surface tension parameters. The upshot of this is 
that our estimates are sufficiently strong to allow us to pass to the vanishing surface tension 
limit, establishing a sort of continuity between the problems with and without surface tension. 
We state this result now. 

Theorem 2.9. Let N > 3 be an integer. Suppose that [[/?]] < 0 and cr± > 0, and assume that 
(q5,v%>Vo) satisfy the compatibility conditions of Theorem, 12.,91 and the mass condition (12.201) 
in addition to the bound £f N (0) + J-^.y(0) < n, where k is as in the theorem. Further assume 
that (j -j- —} 0 and that 

Qo —>■ Qo in H iN (fl), Uq uq in H iN (12), 

Vo Vo in R 4Ar+1 / 2 (£), and '/o'^*Vo — > 0 in R 4JV (E) (2.33) 

as a± —>• 0. 

Then the following hold. 

(1) The triple {qo,uo,r]o) satisfy the compatibility conditions of Theorem 12.21 with a± = 0. 

(2) Let (q a , u a , r/ 7 ) denote the global solutions to (12.101) produced by Theorem 171,71 from, the 
data (qo,Uo>Vo). Then as a± —» 0, the triple (q a ,u a ,r] a ) converges to (q,u,ij), where 
the latter triple is the unique solution to (12.101) with a± = 0 and initial data {qo,uo,rjo). 
The convergence occurs in any space into which the space of triples ( q , u, rj) obeying 
Q 27 v (°°) < 00 compactly embeds. 

2.4. Strategy and plan of the paper. Note that sufficiently regular solutions of (12.101) obey 
the physical energy-dissipation law 

^ (X + R(ypS>J + 9p ( X3 + (9 ) J + VatmV + °+\l 1 + |V*? 7 | 2 

+ J (J-\! 1 + |V*?7| 2 ^ + /fj|B>| 2 + ^|div^| 2 d, = 0, (2.34) 

where 

R±(z) = z [ ^ ^ ds (2.35) 

Jc± s 

for some choice of c± > 0 such that R± > 0. This energy-dissipation equation is too weak on its 
own to even control the Jacobian of our coordinate transformation, and so we must employ a 
higher-regularity framework for global-in-time analysis. Our framework is based on the energy 
and dissipation functionals obtained from linearizing (12.341) around the equilibrium. 

With these functionals in hand, we develop a nonlinear energy method for studying the 
problem (12.101) . This method is similar to that employed in the incompressible setting |8j, [9, 33], 
though of course the compressibility adds many new complications. Below we summarize some 
of the principal features of our analysis. 

Horizontal energy estimates: The basic strategy in any high order energy method is 
to apply derivatives to the equations and use the basic energy-dissipation structure to get 
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estimates. This is complicated by boundary conditions, which impose restrictions on which 
derivatives may be applied; indeed, the differentiated problem must have the same structure 
as the original in order to produce the same energy-dissipation estimate. In a moving domain 
this is particularly annoying, as the acceptable directions vary not only in space but also in 
time. While it’s possible to get around this difficulty in a moving domain (see for example 
the work of Bae [2]), the more common approach is to change coordinates into a flattened 
(or at the very least unmoving) domain. This is what we have done in (12.101) . In the flattened 
framework, the admissible derivatives are the temporal and horizontal spatial ones because they 
preserve the structure of the boundary conditions. As such, the basic energy estimates can only 
control temporal and horizontal derivatives of the solution. Control of vertical derivatives must 
be acquired from other methods, namely elliptic estimates and intermediate energy estimates 
coming from special structure in the equations. 

The trade-off for flattening the coordinate domain is that the operators become far more 
complicated and fail to be constant-coefficient. We are then faced with a choice: work directly 
with the system (12.101) . or rewrite the problem as a perturbation of the corresponding constant- 
coefficient problem. The benefit of the latter approach is that the differential operators commute 
with all derivatives, making the basic energy estimates cleaner. The downside of this approach 
is that the forcing terms are of higher regularity and may be outside of what can be controlled 
by the energy method. Unfortunately, this is what occurs for (12.101) : estimates of the highest- 
order temporal derivatives do not close in the perturbed form because we cannot control the 
highest temporal derivatives of u on the internal interface. Because of this, we are forced to 
work directly with the “geometric form” (12.101) for some of the energy estimates. Note, though, 
that we can use the “perturbed form” for some of the energy estimates, and indeed, this is more 
convenient for some. 

Energy positivity: Ensuring the positivity of the energy is one of key issues in using an 
energy method for global-in-time analysis. In the case of a possible Rayleigh-Taylor instability, 
when [p] > 0, the energy appearing in our analysis does not have a definite sign a priori. It 
is only by imposing a condition on cr_ that we can restore this positivity, and this actually 
reveals the stability criterion: it is exactly the threshold at which the energy can be guaranteed 
to be definite. Actually proving this positivity amounts to using the sharp Poincare inequality 
on T 2 . This in turn requires that we control the average of ??. In the incompressible problem 
this is trivial since the average is constant in time, but in the compressible case, we must use a 
more delicate argument, employing the conservation of mass to control the average with other 
controlled terms. 

Intermediate energy estimates: To go from control of horizontal derivatives to control 
of full derivatives we will employ elliptic estimates. Unfortunately, the horizontal energy and 
dissipation do not control everything necessary for such an elliptic estimate, so we are forced to 
seek intermediate energy estimates. These come from the special dissipative structure that arises 
by taking an appropriate linear combination of the density equation and the third component of 
the velocity equation in (12.101) . To the best of our knowledge, this structure was first exploited 
by Matsumura and Nishida [18] • Using this trick, we get an energy-dissipation equation for 
vertical derivatives of the density perturbation and derive estimates. When combined with 
the basic horizontal estimates and elliptic regularity estimates for a certain one-phase Stokes 
problem, these yield control of intermediate energy and dissipation functionals that control 
nearly all the desired derivatives. 

Comparison estimate: We gain control of all desired derivatives by way of comparison 
estimates. In these we use elliptic regularity estimates to show that the intermediate energy 
and dissipation control the “full” energy and dissipation, up to some ultimately harmless error 
terms. The comparison estimate for the energy follows by applying the elliptic regularity theory 
for a two-phase Lame problem. The comparison estimate for the dissipation is particularly tricky 
because the horizontal dissipation itself provides no control of the free surface functions. Indeed, 
the hardest terms to control are ||g|| 0 and ||^|| 0 . For them, we again use the conservation of 
mass to construct an auxiliary function w that allows us to derive the control of ||g|| 0 and ||?7|| 0 
together. The energy positivity mentioned above again plays a crucial role here. 












12 


JUHI JANG, IAN TICE, AND YANJIN WANG 


Decay estimates and the two-tier scheme: The basic goal in a nonlinear energy method 
is to derive an energy-dissipation equation of the form 

^ -£ + V < £ e T> for some 9 > 0. (2.36) 

In a small-energy regime (verified by the local theory and choice of small data), the right-hand 
term can be absorbed onto the left, resulting in an equation of the form 

(2.37) 

In the event that the dissipation is coercive over the energy space, i.e. V > C£, we may 
immediately deduce exponential decay of £. When we consider (12.101) with surface tension and 
( 7 _ > a c in the case [p] >0, we can essentially prove an estimate of the form (12.361) and the 
coercivity of the dissipation, and we can deduce exponential decay of solutions. 

However, when we neglect surface tension, we run into two major difficulties. The first is 
that we cannot derive an estimate of the form (12.361) because the regularity demands for p in 
the nonlinear forcing terms are half a derivative beyond the control of the energy or dissipation. 
This would be disastrous if it were not for auxiliary estimate for r/ arising from the kinematic 
transport equation. This estimate provides enough regularity to guarantee that the nonlinearity 
is well-defined, but the estimate can only say that this term grows no faster than linearly in time. 
This prevents us from using a small-energy argument to absorb the nonlinear term as above. 
The second difficulty is that the dissipation is not coercive over the energy space. Indeed, we 
find that there is a half-derivative gap between the dissipation and energy for p. This prevents 
us from deducing exponential decay of solutions, even if we manage to derive ([2.37D . 

Our solution to this problem is to implement the two-tier energy method devised by Guo 
and Tice in the analysis of the one-phase incompressible problem The idea is to employ 

two tiers of energies and dissipation, £i , T>i , £h, and T>h , where l stands for low regularity and 
h stands for high regularity. We then aim to prove estimates of the form 

j t S h + V h < (£ h ) e V h + m) x l 2 V h and j£ x + V x < {£ h ) e V u (2.38) 

where the t term results from the transport estimate. If we know that £i decays at a sufficiently 
fast polynomial rate and that £^ is bounded and small, we can deduce from these that 

r°° d 1 

sup£ h + V h < £h(0) and —£i + -T>i < 0. (2.39) 

t Jo at 2 

We still don’t have that T>i > C£i, but we can use an interpolation argument to bound £i < 
(T/j) 1 _ 7 ("D ;) 7 for some 7 € (0,1). Plugging this in leads to an algebraic decay estimate for £1 with 
the rate determined by 7. If the high and low regularity scales are appropriately chosen, this 
scheme of a priori estimates closes and leads to a method of producing global-in-time solutions 
that decay to equilibrium. An interesting feature of this is that the existence of global solutions 
is predicated on the decay of the solutions. 

The rest of the paper is organized as follows. In the remainder of this section we record some 
definitions and bits of terminology. Section [3] contains some preliminary estimates and energy- 
dissipation identities that are useful throughout the paper. Section |4] produces the estimates 
for temporal and horizontal derivatives. Section [5] produces energy-dissipation estimates for 
vertical derivatives of the density perturbation. Section [ 6 ] combines the horizontal and temporal 
estimates with the density perturbation estimates to produce estimates for the intermediate 
energy and dissipation. Section [ 7 ] shows that the intermediate energy and dissipation control 
the full energy and dissipation up to error terms. Section [ 8 ] records the finished a priori estimates. 
Section [9] contains the proofs of Theorems 12.31 and 12.91 

2.5. Definitions and terminology. We now mention some of the definitions, bits of notation, 
and conventions that we will use throughout the paper. 

Einstein summation and constants: We will employ the Einstein convention of summing 
over repeated indices for vector and tensor operations. Throughout the paper C > 0 will denote 
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a generic constant that can depend N, 12, or any of the parameters of the problem except for 
surface tension. We refer to such constants as “universal.” They are allowed to change from one 
inequality to the next. We will employ the notation a < b to mean that a < Cb for a universal 
constant C > 0. 

To track the appearance of constants depending on the surface tension coefficients, we employ 
the notation 


K{<j, [pD 


c(cr) if Ip} > o 
B(a) if [pi < 0, 


(2.40) 


where C(a) > 0 is a constant depending on a (and the other parameters of the problem) and 
B(a) > 0 denotes a constant (also depending on the other parameters) satisfying 


liminf B(a) € (0,oo) (2-41) 

(fT-)- ,<T— )—^0 


and 

limsup - r ^ ^ —t-t < oo for some integer j > 0. (2.42) 

( ff+ , CT _)4.oomax{a+,a_p 

Notice in particular that linear combinations of powers of constants of the form K(a, [p]) 
remain of the form K(cr, [[p]) so that, for example, we may write “equalities” of the form 
1 + K(a, [p]) 5 + y/K(a,lp\) = K(a, [p]). 

Norms: We write H k (Q±) with k > 0 and and H S (E±) with s € R for the usual Sobolev 
spaces. We will typically write H° = L 2 . If we write / € H k (tl), the understanding is that / 
represents the pair f± defined on f2± respectively, and that f± € H k (Q±). We employ the same 
convention on E±. We will refer to the space o27(12) defined as follows: 

oi7(f2) = {v £ H 1 (0) | [u| = 0 on S_ and r_ = 0 on Ej,}. (2.43) 

To avoid notational clutter, we will avoid writing H k (Q) or H k (Y,) in our norms and typically 
write only ||-|| fc , which we actually use to refer to sums 

\\f\\k = ll/+ll// fe (Q + ) + ||/-||#fc(f2_) o r \\f\\k = ll/+llH fe (E + ) + ll/-llH fe (s_) • (2.44) 

Since we will do this for functions defined on both 11 and E, this presents some ambiguity. We 
avoid this by adopting two conventions. First, we assume that functions have natural spaces on 
which they “live.” For example, the functions u, p, q, and fj live on fl, while p lives on E. As 
we proceed in our analysis, we will introduce various auxiliary functions; the spaces they live on 
will always be clear from the context. Second, whenever the norm of a function is computed on 
a space different from the one in which it lives, we will explicitly write the space. This typically 
arises when computing norms of traces onto E of functions that live on 12. 

Occasionally we will need to refer to the product of a norm of rj and a constant that depends 
on ±. To denote this we will write 

7 Ml = 7 + ||» 7 +||Hfc(E + ) + 7 - \\y -\\■ ( 2 - 45 ) 

Derivatives 

We write N = {0,1,2,...} for the collection of non-negative integers. When using space-time 
differential multi-indices, we will write Nf 1+m = {a = (ao,ai, •.. ,a m )} to emphasize that the 
0—index term is related to temporal derivatives. For just spatial derivatives we write N m . For 
a £ N 1 + m we write d a = d^d^ 1 ■ ■ ■ d^ 1 . We define the parabolic counting of such multi-indices 
by writing |a| = 2ao + aq + • • • + a m . We will write V*/ for the horizontal gradient of /, i.e. 
V*/ = d\fei + <92/^2, while V/ will denote the usual full gradient. 

For a given norm ||-|| and an integer k > 0, we introduce the following notation for sums of 
spatial derivatives: 


v*7 


2 


E ll^/ll 2 “d 

ae N 2 

|a:|<fc 



E lld“/ll 2 ■ 

aeN 3 

|a|<fc 


(2.46) 
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The convention we adopt in this notation is that V* refers to only “horizontal” spatial deriva¬ 
tives, while V refers to full spatial derivatives. For space-time derivatives we add bars to our 
notation: 


v?/ 


E ii^n' 


and 


V fc / 




ii3°7ir 


(2.47) 


aGN 1+2 aGN 1+3 

|a|<fc |a|<A; 


We allow for composition of derivatives in this counting scheme in a natural way; for example, 
we write 


v*v*/ 


v*v*/ 


E 


'V./ll 2 = 


E 


ir/ir 


(2.48) 


aeN 2 aeN 2 

|a|<fc l<|a|<fc+l 


3. Preliminaries 

In this section we record some preliminary results that are useful throughout the paper. 
We rewrite the problem (12.101) as a perturbation around the equilibrium, and then we present 
various energy equalities. We conclude with the estimate of certain nonlinear terms. 


3.1. Perturbed formulation around the steady-state — geometric form. We will now 
rephrase the PDE ([2.101) in a perturbation formulation around the steady state solution (p, 0,0). 
We recall our special density perturbation defined by (12.141) . In order to deal with the pressure 
term P{p) = P(p + q + d 3 p 8 ) we introduce the Taylor expansion, by (11,111) . 

P(-p + q + d 3 p0) = P(p) + P\p)(q + d 3 p 8 ) + 77 = P(p) + P'{p)q - gpd + 77, (3.1) 

where the remainder term is defined by the integral form 

rp+q+d3p0 

77 = / (p + q + d 3 pd - z)P"(z) dz. (3.2) 

Jp 

Then one can see that the perturbation q is defined in its special way enjoying the following 
advantage that in 17, 

AijdjP(p) + gp5 i3 = AijdjP(p) + gpAijdj @ 3 

= Aijdj(P(p) + P'{p)q - gpO + 77) + g(p + q + d 3 p6)Aijdj(x 3 + 6 ) 

= Aijdj(P'(p)q ) - gAijdj(pd ) + AijdjTZ + gpAijdjd + g(q + d 3 p8)A i3 + g(q + d 3 pe)A ij d j e 
= Aijdj(P'(p)q) + AijdjTZ + gqA i3 + g{q + d 3 p9)Aijdj6 
= pAijdj(h'(p)q ) + AijdjTZ + g(q + d 3 pd)Aijdj8, 

(3.3) 

where we have used (11.111) and (11.131) . Recalling also ([1.261) . (12.81) and (12.51) . we have 

- 9P+0 = ~Pi9V+ on s +> and 1~9P0 1 = - M 9V- on E_. (3.4) 

By perturbing the density as above, we are led from (12.101) to the following system: 


’ d t q + div^(pu) = F 1 

(p + q + d 3 pd)d t u + pS7 a ( h'{p)q) - div^S^u = F 2 

dtp = u ■ M 

< ( P'(p)ql - Sa(u))M = pigp+N - cr + A* 77 + Af + Ff 
\P'{p)qI - §4(14)]] M- = [p] grj-JV + o\_ 

M = ° 

K U- = 0 


in 17 
in 17 
on E 
on E + 
on E_ 
on E_ 
on E b , 


(3.5) 


F 1 = d\pK 8 d t 8 + Kd t 8 d 3 q - div A ((q + d 3 pd)u), 

F 2 = -(p + q + d 3 p 8 )(-Kd t 8 d 3 u + u ■ V A u) - V A U - g(q + d 3 p8)\7 A 8 , 
Fl = -TZN - cr _j_ div*(((l + jV^+l 2 )- 1 / 2 - l)V*r/+)A f, 


where 


(3.6) 

(3.7) 

(3.8) 
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and 

- Fl = - IK} Af + <r_ div*(((1 + IV^.I 2 )- 1 / 2 - l)V,p_)M. (3.9) 

We will employ the form of the equations (13.51) primarily for estimating the temporal deriva¬ 
tives of the solutions. Applying the temporal differential operator d{ for 7 = 9...., 2N to (13.51) . 
we find that 


' d t (diq) + dw A (pdiu) = F 1 ’! 

in 11 


{p + q + d 3 pd)d t (d{u) + pV A (h'(p)d/q) - div^S. A(d{u) = F 2 ^ 

in hi 


dt(d{p) = d{u ■ Af + F 4j 

on E 



P'(p)d{ql - S A (d/u))Af = pxgdlp+N - a+A^d/p^Af + 

P'(p)d{ql — S> A (d{u) Af = [pj gd{p-Af + cr_ A*( cP t p-)Af — F^ 
d-j-u =0 

on E + 

on E_ 

on E_ 

(3.10) 

dlu- = 0 

L 

on E fe , 



where 

F 1 * = djF 1 - Y CffiA lk d k (pdi~ £ Ul ), (3.11) 

o <e<j 

F i ' 3 = d t F i + C J {^ikdkidfAimdl^dmUi) + pdfAi k di~ i d k (Ai m d m Ui) 
o <t<j 

+ {p/3 + jj!)Ai k d k (dfAimdi d m ui) + (p/3 + p')d/Ai k di d k (Aimdm,ui) (3-12) 
- pdfAi k d k (h'(p)di~ e q ) - df(q + d 3 p9)d t (d{~ £ u )}, i = 1,2,3, 

F?£ = %!%++ Y C'j' {M+^(AA z ^ fc )^- £ a fc ix z + 

o <l<j 

+(p+ - 2n + /3)dt(AT i Ai k )d 3 t ~ l d k ui + d^M i di~ l {p 1 gr] + - P'(p)q - (T + A*7y + )| , i = 1,2,3, 

(3.13) 


_/? 3J _ _ 


d t F i-+ Y C J {dtWiAi^di ^ lpd k u t J + dfiAfiAi^dj 1 {pd k Ui\ 


o <e<j 


+df(A fiAi k )di 1 {(g, 1 - 2/x/3 )d k uij + d l t Afid{ \{pj gp. - \P'(p)q\ + cr_A*? ? _)} , 

(3.14) 

for i = 1,2, 3, and 

F 4j = Y C]d[N ■ dP t ~ l u. (3.15) 

o <£<j 

Note that the equations (13.101) for (djq,d/u,cP t p) for j = (),..., 2N have the same structure 
in their left-hand sides. We present their energy identities more general way as follows: 


Proposition 3.1. Suppose that ( q,u,p ) solve (12.101) and that A,J, and Af are determined 
through p as in (12.61) . (12.71) . and (12.121) . Further suppose that ( Q,v ,() solve 

in fl 
in fl 
on E 

on E_|_ (3.16) 

on E_ 
on E_ 
on T,k. 


' d t Q + di v A (pv) = 5 1 

{p + q + d 3 p0)d t v + pV A (h'(p)Q) - div^ S A (v) = $ 2 
d t ( = v-Af + s 4 

(P'(p)QI - S A (v))Af = pig(+Af ~ cr+A^+Af + 
lP'(p)QI - S,i(u)] M = M 9(-M + v-AA-AT - g 3 _ 

H = o 

V- = 0 


Then 
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1 d_ 

2 dt 


(p + Q + d 3 p9)J\v\ + h'(p) J\Q\" + / pip|C+l +cr + |V*C+|‘ 
In J e+ 


+ - MiHC-l 2 + cr- |l 2 ^ + | J\^A V \ 2 +p'J\divAv \ 2 

= \ I d t (J(p + q + d 3 p 8 ))\v \ 2 + h'(p)d t J\Q\ 2 + f J(ti(j>)Qt+v$ 2 ) 

1 Jn J n 

+ [ -v-t+ f PigC+t + - [ MaC-t-- f crA^a 4 . (3.17) 

i/ S i/ S-|- i/ S i/ s 

Proof. Taking the dot product of the second equation of (13.1611 with Jv and then integrating 
by parts over the domain Q, using the Dirichlet boundary conditions of v, we obtain 

^4 [ {p + Q + d 3 pO)J\v \ 2 + j ^j\O 0 A v \ 2 + p'J\div A v \ 2 

= 77 / d t {J{p + q + d 3 pe))\v \ 2 + f Jh'(p)Qdw A (pv) + f Jv -$ 2 (3.18) 

1 Jn Jn Jn 

- [ ( P\p)QI - S A (v))U • v + [ {P\p)QI - S A (v)j M • v. 

JY. + 

Using the first equation of (13.1611 . we have 

[ Jh'(p)Qdw A (pv)= [ Jh\p)Q(-d t (Q)+t) 

Jn Jn 

= ~\j t f h'(p)J\Q \ 2 + ± j h'(p)d t J\Q \ 2 + f Jh'{p)Qt- 

Using the fourth equation of (13.1611 . we have 

~[ (P\P)Q I ~ §a{v))N ■ v = - [ (piffC+A7-cr + A*C + A7 +5^) • u. 

Jn+ J e + 

The third equation of (|3.16l) further implies 


/ ( Pig(+ - o-+A*C+)A7 • v 
J E + 

= - / (pipC+ - c r +A*C+)(^C+ - 5+) 

Jt.+ 


ld_ 
2 dt 


PiSlC+l +0+|V*C+l + / (PiS<+ - 0+A*C+)3+. 


Hence 


- [ (P\p)QI-§ a {v))M-v 

Jt. + 

= P^lC+l 2 + a +\^*(+\ 2 + ~ v ■ 5+ + (pigC+ - o+A*C+)3+- 

Similarly, we use the fifth and third equations of (13.161) to have 

/ iP'(p)QI-§ A (v)}N-v 


ld_ 
2 dt 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


- M g IC+I 2 + cr- |V*C-| 2 + -v-d- + (- [p] g(- - cr_A*C-)3_. 


(3.23) 


Consequently, plugging (|3.19|) . (13.22|1 and (|3.23l) into (13.181) . we obtain (13.171) . 


□ 
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3.2. Perturbed formulation around the steady-state — linear form. It turns out to be 
convenient to write the system (12.101) in a linear form to derive the subsequent estimates. The 
reason for this is that the operators become constant-coefficient, wl 
elliptic regularity. We rewrite the PDE (13.51) for (q, u, iq) as 

'd t q + div(pu) = G 1 
pdtu + pS7 (h'(p)q) — div§(u) = G 2 
dtp = u 3 + G 4 


lP'(p)qI-S(u)j e 3 = 

M = o 

it- = 0 

where we have written the function G 1 = G 1,1 + G 1,2 for 

G 1 ’ 1 = Kd t 9d 3 q - uiAi k d k q , 


l is 

more 

convenient for 

in 

n 


in 

n 


on 

E 


on 


(3.24) 

on 

£_ 


on 

£_ 


on 




(3.25) 

(3.26) 


(3.27) 


G 1 ’ 2 = d 2 pK0d t 9 - qAikd k ui - Ai k d k (d 3 p9ui ) - (A ik - Si k )d k (pui), 
the vector G 2 for 

Gj = - (q + d 3 pB)d t Ui + (p + q + d 3 p0)(Kd t 9d 3 Ui - uiA lk d k Ui) 

T P j Ai k d k Ai m d rr iUi -{- p(A[ k Ai m ^ik^im^)d km Ui 

T (p/3 T p )Ai k d k Ai m d rn ui + (p/3 + p ) (Ai k Ai m ^ik^im)d k mUi 
~ P(Al ~ Su)di(ti(p)q) - Aud{R - g(q + d 3 p9)Audi9, i = 1,2,3, 

the vector G+ = G^ 1 + &+G+ 2 for 

Gf >+ = p + (Audiu k + A k idiUi)(N k — 5 k3 ) + P+(Au — 5u)diu 3 + p(A 3 i — 5 3 i)diUi 

+ (p' + - 2p + /3)Ai k d k ui(Ni - 5 i3 ) + (p' + - 2p + /3)(Ai k ~ 5i k )d k ui5 i3 + p 1 gij + (M i - < 5 ; 3 ) 
-nMi + P'(p)q(5 i3 -Ni) 

(3.28) 

and 

Gj 2 = -A* v+ (M t - 5 a ) - div*(((l + IV^+I 2 )- 1 / 2 - 1)V*7 1+ )Mi (3.29) 

for i = 1,2, 3, and the vector G :i _ = G^) 1 + <t_G^ 2 for 

-G- _ = (An IpdiUkj + A ki \pdiUi\)(N k - S k3 ) + (A a - S a ) \pdiu 3 } - (A 3i - 5 3l ) [[ pd t Ui} 

+ Aik l(p' - 2p/3)d k uij (J\fi — d i3 ) + (Aik ~ <5ik ) {(p ~ 2p/3)d k uij S i3 
+ [Pi 9V-(M’i - S i3 ) - pi} Mi + {P'(p)qj (S i3 - Mi), 


and 

G 3 ! 2 = Mp-(Mi - S i3 ) + div*(((l + |V*t/_| 2 ) -1 / 2 - 1 )V*ri-)Mi 

for i = 1,2, 3, and the function G 4 for 

G 4 = -uidip - u 2 d 2 p. 

In all of these, H is as defined in (1 3. 2D . 

We now present the energy identities related to (13.241) . 


(3.30) 

(3.31) 


(3.32) 
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Proposition 3.2. Suppose that (Q,v,() solve 

' d t Q + div(/w) = © 4 
pd t v + pV ( h'(p)Q ) — div§(u) = © 2 

dtC = V 3 + © 4 

< ( P'(p)QI - §(u))e 3 = (pig(+ - cr+A*C+)e 3 + ©+ 

[- P'{P)QI ~ %)! e 3 = ([pl g(- + cr_A*C-)e 3 - © 3 

M = 0 

K V- = 0 

Then 


in ft 
in ft 
on £ 
on £+ 
on li¬ 
on li¬ 
on £;,• 


(3.33) 


\^ t h'{p) \Q \ 2 + p\v \ 2 + pig |C+| 2 + cr+ |V*C+| 2 


+ / ~ M 9 IC-I 2 + <?- |v*C-l 2 ) + 


'E_ 


T i™o |2 


l Z i / I A- L 

v\ + fi |dm;| 


= [ ti(p)Q&- +vd a ® 2 + [ -v & 3 

Jn J e 

+ [ Pig(+®\+ f -[pIK-©-- [ dA*C© 4 - (3.34) 

J S_i_ 1 S J S 


Proof. The proof proceeds along the same lines as that of Proposition 13.11 we multiply the 
second equation in (I3.33P by v, integrate over and integrate by parts and the sum 
resulting equations. After employing the other equations in (|3.33l) as in Proposition 13.11 we 
arrive at (13*341) . □ 


3.3. Estimates of the nonlinearities. We assume throughout this subsection that the solu¬ 
tions obey the estimate Qf N (T) < 5, where 6 € (0,1) is given in Lemma I A. 31 

We first present the estimates of the nonlinear terms G 1 (defined by (I3.25l) - (|3.32p f at the 2 N 
level. 


Lemma 3.3. It holds that 


V 4 


2 g 1 ||^ +11 v 4 


G l|o + ll V * 


2 g 3 ii ; /2 + iiv: 


4 g 4 


11/2 


\P\)^ 2 N^ 2 N + £% + 2 ^ 2 N, (3.35) 


and 


|y47V—1 G 1,1||2 + || V 


4Af - 2 d t G 1 ’ 1 || 2 + ||V 4 JV G 1 ’ z ||q + ||V 4JV - 4 G : 


74Ar^l,2|| 2 


4AT-1/-<2|| 2 
0 


+ || V 4 Ar - 1 G 3 ||^ /2 + || V ^- 1 G 4 ||^ /2 


+ II Vi N ~ 2 d t G 


+ a 2 ||V 4A, G 4 


li/2 


1/2 

2 N^ 2 N + £n+ 2 ^ 2 N- (3.36) 


Proof. We first prove the estimates in (13.361) . Note that all terms in the definitions of G l are at 
least quadratic. We apply these space-time differential operators to G l and then expand using 
the Leibniz rule; each term in the resulting sum is also at least quadratic. We then estimate 
one term in H k {k = 0 or 1/2 depending on G 1 ) and the other terms in H m for m depending on 
k, using trace theory, and Lemmas IA.lllA.2l along with the definitions of Sf N and P% N ( dU5]> 
and (|2.16l) . respectively), and K(a,\p\). With three exceptions, we can estimate the desired 
norms of all the resulting terms by £ 2Ar D 2Ar . 

The first exceptional terms are ones involving either V 4Ar+1 ?7 in kl or V 4Ar ?? on £ when 
estimating ||V 4iV G 1,2 || 2 , ||V 4JV_1 G 2 || 2 , || V 4Ar ~ L G 3,1 || 2 / 2 and || V 4A, ~ 1 G' 4 || 2 / 2 . But we can argue 
as Theorem 3.2 of jS] to estimate ||f?|| 2 jv+i ~ ll^llLv+i^ ~ J~ 2 N to bound these terms by 
£n+2 ^~2N- 
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The second exceptional terms result from the gG 3,2 term when estimating || 1 G 3 ||^ 9 
and ||Vf Ar ~ 3 9 t G 3 ||' / ,,- ) . The highest derivative appearing in gG 3 ' 2 is crV 2 ?/. We make use of this 


a factor to estimate 

7 4AT-1v72 m2 




V*»71| 1/2 < O* MIn+3/2 = miu { L(T 2} min {!w 2 } IMlLv+3/2 

nN<K{aMW 2 N, (3.37) 


< 


and, similarly, 


hvi N ~ 3 vld tV \ 

Using these, we may estimate 


2 

1/2 


min{l, a 2 } 


V 2 IMlLv-1/2 ~ K (.°> M) V 2N- 


V 


4V—1 


vG 3 ’% 2 <K(aM)£%V° 2N . 


(3.38) 

(3.39) 


The last exceptional term is the term a 2 || V* Ar+1 ry||^ 2 that appear when we estimate the term 

g 2 ||V 4iV G 4 ||y 2 . Again, the factor er 2 leads to the estimate a 2 ||V 4JV+1 ? 7 ||y 2 ~ 7f(cr, lp})'D% N - 
Hence, in light of the above analysis, we may deduce (|3.36l) . The proof of (I3.35|) proceeds 
similarly. We remark that the term £^ + 2 F 2 n in (13.351) does not appear in Theorem 3.2 of [5]; 


it appears here because we want to control V 


4V-l£<4|| 2 


' 1 / 2 ' 


□ 


We then present the estimates of G l at the iV + 2 level. 
Lemma 3.4. It holds that 


y2(V+2)-2 G l 

2 

+ 

y2(V+2)-2 G 2 

2 

+ 

v 2(V+2)-2 G 3 

2 

+ 

v 2(V+2)-l G 4 


1 


0 


1/2 



< 


2 

1/2 

0 c-0 


KM p])£ 2 Vv+ 2 , (3.40) 


and 


t^2(N+2)q! 

2 

+ 

■^j2(N + 2 )—l q 2 

2 

+ 

V*(^ V_ I~2) — 1(//3 

2 

+ 

V* (jv+2) G 4 


0 


0 


1/2 

< K(a, I 


2 

1/2 


A>,[p])£ 2 V^ +2 . (3.41) 


Proof. The estimates (13.401) (13.4111 follow from arguments similar to those used in the proof of 
Proposition ^. 31 In this case they are easier because when we estimate the terms appearing from 
the Leibniz rule expansions we do not have exceptional terms as in the proof of Proposition 13.31 
Indeed, we may write each term in the form XY, where X involves fewer temporal derivatives 
than Y ; then we simply bound the various norms of Y by K(a, [p])£ 2 iv an d bound the various 
norms of X by 8% +2 for (13.401) or V % +2 for (13.411) . Note here that we use S 2N and absorb all 
appearances of a into K(a, [p]]) rather than employ Sf N - HI 

Next we present some variants of these estimate involving integrals. First we consider prod¬ 
ucts with derivatives of G A . 


Lemma 3.5. Let a € N 2 so that lal = 41V. Then 


L 


d a r]d a G 4 


^2N^2N + \/'^ > 2N^N- t-2^ r 2V 


(3.42) 


and 


L 


aA*d a r]d a G 4 


K(g, + K{a, M)^V° n 8% +2 F 2N . (3.43) 
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Proof. We first use the Leibniz rule to expand that 

- 8 a G 4 = d a {V,r]-u) = V*3V« + XI C^*d a ~ P vd p 

0 </3<a 

For the second part, we estimate for |/?| > 1, similarly as in Lemma 13.31 


u 


V*(9 a-|8 ?7 • d !3 u 


1/2 


< £‘2n'E > 2N + £n+2-£2N- 


(3.44) 


(3.45) 


Hence, we have 

<9“t/V*<9“ _/3 ?7 • d p u < ||<9 a 77||_ 1/2 • d p u < y / £>X \J 8 2N V 2N + &N+2?™ 

(3.46) 


and 


J aA*d a rjV*d a -Pri -O^u^a \\A,d a r ]\\_ 1/2 ||v*<9 a -^ • 8 ? 

mta ( bl }\ /min{1 - <r2)ll ’ ,ll “+ 3 / 2 


X7*d a ~Pri ■ d 0 u (3.47) 


1/2 


< 

rsj 


K(<T, m)V^\l^2N V 2N + ^ + 2^V- 

For the first term, we use integration by parts to see that, by Lemma lA.21 

f d a r)’V*d a r) ■ u = ^ / V* \8 a p\ 2 • u = — ^ f \d a rj\ 2 div* u 
J s 4 z ,/e 

11^11-1/2 ll^“ ? ?lll/2 ll div * U \\h2(Z) < ^ V 2N :F ‘2n£n+2 


(3.48) 


< 

rs_/ 


and 


cr A*8 a rjS7 *8 a p ■ u = — X / &8id a r)8i(V*d a ri ■ u ) 


2— 1 


2 a 


[ V* |V*<9“?y| 2 • u - X [ vdid a r]('V*d a ri ■ di.u) 
4s i=1 4s 

X f |V*0“r/| 2 div* u — X / crdid a p(\/^ 8 a p ■ diU ) 


(3.49) 


2—1 


< a ||V*d°y|| 1/2 ||d a V*7/||_ 1/2 ||V,«|| h2(s) < tf(a, [p])^^ 2 ^ +2 . 
The estimates ()3.42|1 (13.43|) then follow by summing the estimates. 

Next we consider products with derivatives of G 1,1 . 

Lemma 3.6. Let a € N 3 so that |a| = 4 N. Let f(p ) denote either 1 or h'{p). Then 

Jj(p)d a qd a G^ < /^XyX 0 A + ^ + 2 ^v- 
Proof. We hrst use the Leibniz rule to expand that 


□ 


(3.50) 


S^G 1 ’ 1 = LT&tfcl 


' 3 d a q-UjAj k d k d a q+ X (d 13 (Kd t 9)8 3 d a 13 q - (ujAjk)d k d a 13 q^j . 


0</3<a 

For the second term, we estimate for |/3| > 1, similarly as Lemma 13.31 

d^KdtOW^q - d^ 3 (ujA.jk)d k d a ~^q 2 < L 2 %P 2 V + 4 +2 ^2iv. 


(3.51) 

(3.52) 
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Hence, we have 


J f(p)d a q (d fi {Kd t 9)d 3 d Q ~ 0 q - d^ Uj A jk )d k d a ^q 
dP(Kd t 0)d 3 d a -Pq - d' 3 ( U] A Jk )d k d a -' 3 q 


< Iloilo 


~ \h^2N\ + £ V+2^2V- 


~ V 2N V 2 N^ 2 N ^ ^N+ 2 “ 

For the first term, we observe that 


f(p)d a q(Kd t 9 d 3 d a q - UjA jk d k d a q ) 


/(p)iL<9 t 0d 3 |d a <2f - 


/n 


We then use the integration by parts to have 

f(p)Kd t 9 d 3 \d a q \ 2 - f(p)ujA jk d k \d a q \ 2 


' n_i_ 


f(pi)(Kd t e - UjA j 3 )\d a q \ 2 - f f(p + )(Kd t 9 - UjA j 3 )\d a q \ 2 


{d 3 {f{p)Kd t 9) - d k {f{p)u j A jk )) |<9 Q <?| 2 


(■ d 3 (f{p)Kd t 9) - dk(f{p)ujAjk)) \d a q\ 2 


< 


|| d 3 (f(p)Kd t 9) - d k {ujAjk) 


v|| 2 < /rO -r)0 
9llo ~ \J c '2N U 2N 


and 


f f{p)Kd t 9 d 3 \d a q \ 2 - f{p)u j A jk d k \d a q\ 

Jn- 

f(p~){Kd t 9 — UjAj 3 )\d a q -\ 2 — [ {d 3 (f (p)Kd t 9 ) - d k {ujA jk ))\d a q\ 

e_ Jn 


(d 3 (f(p)Kd t 9 ) - dkiujAjk)) \d a q \ 2 


< 


(3.53) 


(3.54) 


(3.55) 


(3.56) 


\\d 3 (f(p)Kd t 9) - d k ( Uj A jk )\\ Loo \\d a q\\ 2 0 < ^Jsf N V%. 

Here we have used the facts that UjAj 3 = KujAfj = K'dtp = Kdt9 on X, and 0 = 0 and U- = 0 
on Xfe. 

Then the estimate (13.5011 follows from (13.5311 and (13.5511 . □ 

Next we consider a similar estimate involving weights and derivatives of G 1,1 . 

Lemma 3.7. Let a € N 1+3 with |a| < 41V and ao < 2 N — 1. Then 
4/x/3 + p n 


1 + j d a (h'(p)q)d a (h'(p)G 1 ’ 1 ) < K(a, m)y/V^y/S^2N + ^N+2^ N . 

(3.57) 

Proof. We split to cases. For the case 1 < ao < 21V — 1 or the case that «o = 0 and |a| < 41V — 1, 
we use the estimates (13.3611 of Lemma 13.31 to bound 

(i + 4 y| i) + / ) trih’mtr (fc'WG 1 ' 1 ) 

< ||V 4 M| 0 (||V 47V_2 G 1 ’ 1 || 1 + ||v 4JV_3 9 t G 1 ’ 1 || 1 ) 

K(t M)y/vy\jK{v, IpW$ n VZn + S° n+2 ^2 n- (3.58) 


< 
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For the remaining case in which ao = 0 and |a| = 4 N, we rewrite 

h'(p)G 1 ’ 1 = Kd t 9d 3 (h' (p)q ) - UjA jk d k {ti(p)q) - Kd t dd 3 (ti(p))q + UjA j 3 d 3 (ti(p))q. (3.59) 

Write Z\ for the sum of the first two terms on the right of (13.591) and Z 2 for the second two. 
We may argue as in the proof of Lemma 13.31 to estimate 

4/u/3 + //' 


in 


1 + 


d a (h'(p)q)d a (Z 2 ) 


h'(p)p 2 

For the Z\ term we must utilize the total derivative structure 
(u< (^\ „\ ns & aaaan,u^\„\ „. a 

12 


\J^' 2 n'^ > 2 N + £n+ 2 ^~ 2 N ■ (3.60) 


d a (ti(p)q)(Kd t 9d 3 d a (ti(p)q) - u,A jk d k d a {ti(p)q) 

I d a {ti(p)q)i 


= KdtOd■ 


iij Aj k 0 k 


\d a {h'{p)q)\ z 


This allows us to argue as in the proof of Lemma 13.61 to bound 
4///3 + // 


1 + 


h’(p)p 2 


d a (ti(p)q)d a (Zi) 


< 


k M ) y/^2N \/^2 n'^ > 2N + ^N+ 2 ^~2N- 


(3.61) 


(3.62) 


Combining the Z\ and Z 2 estimates then gives the desired estimate when ao = 0 and |a| = 41V. 

□ 

We next present the estimates of the nonlinear terms F % ' 3 , defined by (13.111) ([3.151) . at the 
2 N level. 

Lemma 3.8. Let F 1,3 be defined by (13.1 1 D ([3.151) for 0 < j < 2 N. Then 


|i? 1 J|| 2 _i_ ||F 2 ’ J '|| 2 + ||F 3j '|| 2 + \\F 4 ’ j \\ < K( 




rV°o 


(3.63) 


Proof. Note that all terms in the definitions of F 1 ’ 3 are at least quadratic; each term can be 
written in the form XY, where X involves fewer temporal derivatives than Y. We may use the 
usual Sobolev embeddings, trace theory and Lemma lA.il along with the definitions of £ 2N and 
V" N to estimate \\X\\ 2 LOO < £% and || y || 2 < K(a,lpj)V° 2N . Then \\XY\\ 2 0 < \\X\\ 2 LOO || y || 2 < 
K(cr , \P\)£ 2 n'^ > 2 Ni an d the estimate ([3.631) follows by summing. □ 

Next, we present the estimates of F l)3 at the N + 2 level. 

Lemma 3.9. Let F 1,3 be defined by (13.111) (13.151) for 0 < j < N + 2, we have 

Iloilo + Iloilo + ll^'llo + Iloilo £ [Pl)^2°IV ^ +2 


(3.64) 


Proof. The proof proceeds similarly as Lemma 13.81 with the easy replacement of ||y||g < £ 2N 
and ||X||g < Pjv_|_20. As in Lemma [3~T1 we absorb any appearances of a into K(a, [p]) so that 
we can use both £ 2N and V° N+2 . □ 

Finally, we present an estimate of an auxiliary nonlinearity that will be useful later. Define 

$± = (q± + d 3 p±9)d 3 9. (3.65) 

Lemma 3.10. Let be as in ()3.65l) . Then for n = 2N or N + 2 we have 


Ep 

3=0 


< 


£ n mm{£ n ,V n } 


Proof. The proof is straightforward, so we omit it for the sake of brevity. 


(3.66) 

□ 


4. Horizontal energy evolution 

In this section we derive estimates for temporal and horizontal spatial derivatives. We assume 
throughout this section that the solutions obey the estimate Q 2 N (T) < 5, where 5 £ (0,1) is 
given in Lemma lA.31 
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4.1. Energy evolution in geometric form. We now estimate the energy evolution of the 
pure temporal derivatives at the 2 N level. 


Proposition 4.1. For 0 < j < 2 N, we have 


d_ 

dt 


(p + q + d 3 pd)J 


dju 


+ h'{p)J djq 


+ / pig 

jE+ 


djv+ 


+ CT+ 


V*d t J ii + 


+ ~\p\g dtV- + 0 -- +// 


div dju 


K(a, 

(4.1) 


Proof. Applying dj to (12.1011 leads to (13.101) . which is of the same form as (13.161) with v = dju, 
Q = djq, ( = djg, and = F k ’ 3 for k = 1,2, 3. Then Proposition 13. II yields the equation 


~ ( / (P + 9 + d 3 p#)J 


2 dt 


in 

+ 


dju 


+ h'{p)J djq 


-I pig div- + a - V*div- 1 


+ / Plff 

+ /J j 


djv+ 




+ o-+ 


V*<9(r/+ 


+ p' J div .4 Of u 


1 


2 

+ / -$u-F 3 ’ J ' + 


d t (J(p + g + <9 3 /9#)) <9 {u +h'(p)d t J djq + / J(h'(p)djqF 1,3 + dju ■ F 2 ’ 3 ) 


L 


f pigdjq+F ^ 3 - f \p\gdjg-F ^ 3 - j aA*(djr])F 4 ’ 3 . (4.2) 

J j E- i/E 


We now estimate the right hand side of (14.21) for 0 < j < 2N. For the first two terms, we 


may bound as usual ||<9t J|| 

Z/ 00 < and \\d t {J(p + q + d 3 pd))\\ Loo < \f£f N to have 


1 


2 Jn 


d t (J(p + q + d 3 pd )) 


dju 


+ h'(p)d t J djq 


2 N 


dju 


+ 


djq 


(4.3) 


< 


V / ^2JV^ > : 


2V- 


By Lemma 13.81 we may bound the F ] ’ 3 and F 2,3 terms as 


J J{h\p) 8 P t qF^ + dju ■ F 2 ’ 3 ) < \\djq\\ o ||E lj '|| 0 + 


dju 


I f 2 ’ 3 I 


<9L(a,[pl)V^V^^- ( 4 - 4 ) 

For the i 7 " 3 ’- 7 and i 7 " 4 ’- 7 terms, by Lemma [3. 8 1 and trace theory, we have 

I — Q3 U . _|_ I TJ-^-FjP^ ‘ — IF nil nr)j n — / rr A (rtf 


t / piy^'/+-+ t [ -Mgdjv-F^ - I crA *(djr))F 4 

/S JE+ JE 


< 


< 


dju 


||E 3j ]| n + <7 

R-°(E) 11 110 


djv 


\F 4 ’ j \ 


(4.5) 


dju 


+ 


djv ) A (°"> EpH ) V /< ^2iV-^ > 27V ~ ^ (°"> M) \Z^ 2 N \/^ 2 N^ 2 N m 


Finally, we may argue similarly as in Proposition 4.3 of [8], utilizing Lemma lA.31 to get 




Kdj* 


P 


~ L 2 


B°dju 


+ P 


+ p!J div^4 dju 

We may then combine (14.31) (14.61) to deduce (14.1|) from (14.2|) 
We then record a similar result at the IV + 2 level. 
Proposition 4.2. For 0 < j < N + 2, we have 


div dju 


- (4.6) 

□ 


d_ 

dt 


( [ (p + q + d 3 p9)J dju +h\p)J djq + [ p 3 g djp + + cr + V*djr) + 
\Jn Jt.+ 
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+ / ~Ma fyri- 


+ cr_ 


V*<9 {p- 


+ 




/n 




< 


+ a 


div d{ u 


K{aM)V^N V N+ 2- (4-7) 

Proof. The proof proceeds along similar lines as the proof of Proposition 14.11 using instead 
Lemma 13.91 and replacing T>f N by 'D < f J+2 accordingly in the estimates (14.311 and (14.61) . □ 

4.2. Energy evolution in the linear form. To derive the energy evolution of the mixed 
horizontal space-time derivatives of the solution we will use the linear formulation (13.2411 . 

We now estimate the energy evolution of the mixed horizontal space-time derivatives at the 
2N level. 


d_ 

dt 


Proposition 4.3. Let a £ N 1+2 so that ao < 21V — 1 and |a| < 41V. Then 

p\d a u \ 2 + h\p)\d a q \ 2 + [ pig\d a ri + \ 2 + cr + \\7*d a r ] + \ 2 
~l JT.+ 

+ [ ~Mg\d a r]-\ 2 + a-\V^d a r 1 -A + [ ^\n°d a u \ 2 + p'\divd a u \ 2 ( 4 - 8 ) 

je_ ) Jn 2 

K i 0 ': Ip] ) \Z^2n'^ > 2N + (cr, I/l] ) \JtJ > 2N^'N+2^ T ‘ 2 N ■ 

Proof. Since all the three boundaries of Ll are horizontally flat we are free to take the horizontal 
space-time differential operator in the system (13.241) . Applying d a to (|3.24l) . we arrive at an 
equation of the form (13.33[) with v = d a u, Q = d a q, ( = d a p, <5 k = d a G k for k = 1, 2,3. Then 
Proposition 13.21 provides us with the equality 


Id 
2 dt 


(I h'(p)\d a q \ 2 + p\d a u \ 2 + [ pig \d a p +\ 2 + cr + |V*<9"r/ + | 2 

\Jn J s+ 

+ [ -lp}g\d a ri-\ 2 + a_\V*d a r l -\ 2 ') + [ £ |D°«9“u| 2 + y! |div9 Q u 

J Jn 2 

= [ h'(p)d a qd a G l + d a u ■ d a G 2 + / -d a u ■ d a G 3 

Jn J s 


(4.9) 


+ 


[ Pigd a rj + d a G 4 + + [ -[ p}gd a V -d a Gi- f aA*(d a r,)d a G\ 

j s_i_ i/s_ j yi 


We now estimate the right hand side of (14.91) . We first estimate the G 2 ,G 3 ,G 4 terms. We 
assume initially that |a| < 41V — 1. Then by the estimates (13.361) of Lemma 13.31 we have 


/ 

Jn 


d a u • d a G 2 


< \\d a u\\ 0 ||5“G 2 || 0 < ^D^^/K(a, M)£f N Vf N + £% +2 T 2N . 
Similarly, by the estimates (13.361) of Lemma 13.31 along with trace theory, we have 
d Q u ■ d a G 3 < H^ll^ofs) ||5"G 3 || 0 < ||c? a u||! ||9 Q G' 3 || 0 
\/'^2N\j^ (°"> \p\)^2N^2N 4 “ ^N-\-2^~2N 


(4.10) 


L 


(4.11) 


< 

rs_/ 


and 


f Pl gd a r 1+ d a G\+ f -lp}gd a V -d a Gi- [ aA*(d a r,)d a G 4 

J S_|_ </ S_ J s 

< (Iloilo + o- ||5 a r?|| 2 ) ||9 Q G 4 || 0 < K(a, I pi)V^^2N^2N + £^2 n- 


(4.12) 


Now we assume that |a| = 41V. We first estimate the G 2 , G 3 terms. Since op £ 2 N — 1, then 
d a involves at least two spatial derivatives, we may write a = /3 + (a — /3) for some /3 € N 2 with 
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|/3| = 1. We then integrate by parts and use the estimates (I3.36P of Lemma [3731 to have 


/ 

Jn 


d a u • d a G 2 


/ 

Jn 


d a+ 0 u . d a ~ 0 G 2 


QOt+P 


u 


d a ~ p G 2 


< ll^ull, ||vr" 1 G' 2 || 0 < m)£ U 2N^2N + ^N + 2^N. 

Similarly, by using additionally trace theory, we have 

! G 3 | < ||c 

|y 4 IV- 1 G 3 | 


(4.13) 


d a G 3 

= 

[ d a+ 0 u-d a ~ 0 G 3 

< 

d a+0 u 


d a ~ 0 G 3 



Jn 



H- 1 / 2 (S) 



1/2 


< 

rs_/ 


H 1 / 2 (E) 


| 1/2 < ll^llr ||v: 

\/'^ > 2N\j^ ( <T ’ \Pi)^2N^2N + £n+2-F‘ 2N - 


AN —1 q3 


11/2 


(4.14) 


We then estimate the G 4 term, splitting to two cases: ao > 1 and op = 0. If olq > 1> then d a 
involves at least one temporal derivative so that ||<9“r/|| 3 y 2 < \\d^° p \\ 4N _ 2ao+3 /2 < V^ N - This 
together with the estimates (13.3611 of Lemma 13.31 implies 


[ p 1 gd a r 1 + d a G\+ [ -lp}gd a r,-d a Gi- [ a A ,(d a p)d a G 4 
•J E+ i/ E_ i/ S 

< (ll^llo + * ll^ll 3/2 ) II« 9 “ G 4 Hi/2 £ M)V^N^2N V 2N + £ a N + 2^N. 

If cco = 0, we must resort to the special estimates (13.42H (13.4311 of Lemma [331 to have 

[ Pl gd a g + d a G% + [ - [p]| gd a rj-d a G 4 _ - [ aA if d a pd a G 4 
J E+ J E_ 

K(a, M)y/V^^]K(a, M)£% N V 2N + ^N+ 2-^2 V- 


(4.15) 


' s+ 
< 


(4.16) 


We now turn back to estimate the G 1 term, and we recall G 1 = G 1,1 + G 1,2 . For the G 1,2 
part, it follows directly from the estimates (13.361) of Lemma 13.31 that 


/ 

Jn 


h'(p)d a qd a G 1,2 


< 


9llo 


|3“G 1,2 |L < 


V^2N \J^2N^2N + ^N+ 2-^2V- (4-17) 


Now for the G 1,1 term we must split to two cases: ao > 1 and ao = 0. If ao > 1, then by the 
estimates (13.361) of Lemma 13.31 we have 


[ h'(p)d a qd a G U 

Jn 


< 


QWo ll 5 “ Gl,1 Ho £ V^^ 2 N V 2 N + £n + 2 ^ 2 N. (4.18) 


If ao = 0, we must resort to the special estimates (13.5011 with f(p) = h'(p ) of Lemma [3761 to 
have 

h'(p)d Q qd a G hl 


[ 

Jn 


\/'^ > 2N\J^2N^ > 2N + ^N+2^2N- 


(4.19) 


Hence, in light of (14.1011 (14.191) . we deduce from (14.911 that 

(.In P + ^° q ^ + X Pl9 + a+ \ W * dar]+ \ 2 

+ j ~M 9 \d a V-\ 2 + <r_ |V*5“r7_| 2 ^ + J | \m a u \ 2 + p! |divd"u| 


(4.20) 


< 


< 


A>, [p])V^ 


2N V/ ^2N^2N + ^N+2^2N 


K (°"j M ) V^2N^2N + Tv (cr, [p] ) \j'JJ > 2 N^'N+ 2 ^~ w ) 


which is (14.811 . 


□ 


We then record a similar result at the N + 2 level. 
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Proposition 4.4. Let a € N 1+2 so that «o < N + 1 and |o| < 2(N + 2). Then 


d 


dt 


- 7 i ( / p\d a u\ 2 + h'{p) \d a q\ 2 + I p 1 g\d a T ]+ \ z + cr + \S7*d a r] + \ 


|2 


aa' 12 


+ f -Ipjg \ d a r)-\ 2 + (j— |V*d“7 ? _| 2N ) + [ % |lD 0 9 Q Jt| 2 + g! |divd“u| 2 

/ is] 2 


(4.21) 


< 


\p\)y/Z%N V N+2- 


Proof. The proof is the same as that of Proposition 14.31 except that we instead use the estimates 
(13.4111 of Lemma 13.41 □ 

4.3. Energy positivity. We will now verify the key issue that the energy expressions in the 
previous two subsections are positive. 

We first prove the following lemma that provides the value of the critical surface tension value 
a c defined by (12.1911 . 

Lemma 4.5. Suppose that [[/?] > 0, which means that a c . defined by (j2.191) . is non-negative. If 
a- > a c , then the following hold for all C satisfying J t2 C = 0. 

(1) We have the estimate 

(a~ - <r c ) HCll? < o- ||V*C||o - [pj 9 IICIIo • 

(2) If C satisfies 

- cr„A*C - [pl g( = (p on T 2 , 

then we have for r > 2, 


G- - G r 


I r 11 V'' 11 r —2 ’ 


(4.22) 

(4.23) 

(4.24) 


Proof. Note first that 


1 


I |2 • r t —2 t — 2 

\n\ = mm{L 1 ,L 2 } =- rr2 r2l ■ 

n£(L- L z)x(L~ L z) \{o} max{Lf,L^| 


mm 


(4.25) 


Then since f T2 ( = 0, i.e. £(0) = 0, we may use the Parseval theorem to prove Poincare’s 
inequality with the precise constant: 

HV.CIIo = E M 2 |C(n )| 2 

ne(L~ 1 Z)x(L~ 1 Z)\{0} 

1 E IC (")| 2 = 1 


> 


(4.26) 


(4.27) 


ma x{L‘\,L%\ ^ max{L?, Li} 0 

n£(L^Z)x(L^Z)\{ 0} 1 2 

This inequality is clearly sharp. 

From (14.2611 we see that 

II V*Cllo - MS'IICIIo > - [pJpmax{L?,l|}) ||V*Cllo 

= (cr_ -a c ) ||V*CHo- 

Combining (I4.26|) with (14.271) then yields (|4.22|) . 

To prove (|4.24|) . we note that (14.231) implies that the Fourier coefficients of C and <p satisfy 

(4.28) 

(4.29) 

(4.30) 


(<t_ \n\ — [pj g)C,{n) = <p( n ) f° r all n ^ (^i i ^) x (-^2 i ^ 2 ) 
which in particular implies that <p(0) = 0. Again (14.2511 implies that 


cr_ \n \ 2 -{pjg = a _ |n| 2 - 


G r . 


max{L 2 , L^} 


> (cr_ — g c ) |n| 2 > 0. 


Then 


(g_ - G c ) 2 \n \ 4 ((n) < (<T_ |n| 2 - [p]s)((n) = \<p(n)\ z , 


and hence 
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(ct_ - a c 


lr ~ ( a ~ ~ a c) 


E 


(1 + |n| 2 ) r C(n) 


n&{L^ 1 Z)x(L^ 1 Z)\{ 0} 


£ E 

ne(L- 1 Z)x(L~ 1 Z)\{ 0} 


(1 + |n| 2 ) r 


n 


I v(n)| S 


< 


E 


(1 + |n| 2 ) r_2 |<p(n 


ne(L~ 1 Z)x(L- 1 Z)\{0} 

<IMlE- (4-31) 
□ 

We now show the following positivity: 

Proposition 4.6. Let a € N 1+2 so that |a| < 2 n with n = 2 N or n = N + 2. If [p] < 0, then 

i2 


h'(p) |<9“g| 2 + f pig\d a r] + \ 2 + cr_j_ |V*9“r? + | 2 + f - Jp] g |d“p -| 2 + <t_ |V*<9"?7_| 
./s+ ./£_ 

^ll^ + ll^ + allV^r/ll 2 . 


(4.32) 


If [p] > 0 f/ien 

/ ^(p) l^“g| 2 + / Pi5|5“t? + | 2 + ct+ |V*d“r/ + | 2 + [ -lp}g\d a r]-\ 2 + cr_\\7*d a ri-\ 2 
J £1 JE_|_ «/S_ 

> ||<9 a <?||o + min{l,o+} ||d a ? 7 + || 2 + min{l, ct_ - <t c } ||<9 a ?7-||? - y/S^mm{S^, V°} 

(4.33) 

Proof. In the case [[p] < 0 the result (14.321) is trivial, so we assume in what follows that [p]| > 0 
and focus on (14.331) . 

For a with or + 02 / 0, we can easily get (14.331) by (14.221) since in this case £ = d a g- has 
the zero average over E_, and so Lemma 14.51 is applicable. Now for or = 02 = 0, that is, the 
pure temporal derivative case, and the problem is that they do not have the zero average. 

We first deal with the case ccq = 0. We will make use of the conservation of mass: 


This implies 


[ P±J = [ (p± + q± + d 3 p±0)(l + d 3 0) = f p±. 

J H-j- J J £l± 

f q±+ [ d 3 {p±0) + f 4> ± = 0 , 

J £i± J q± J 


where we have denoted = (q± + d 3 p±0)d 3 0 as in (13.651) . On the other hand, we have 


and 


Then we obtain 


/ d3(p+0)=pi r)+-p + g- 

J £1- 1_ «/S_|_ «/S_ 

d 3 (p-0) = p- 

q++ I $+ = -pi / v+ + p 1 ' 

J Q_i_ J E_i_ 


/ 

Jci- 


V- 


I £_ 


V— 


and 


/ ,-W 

J il - J E. 


g~, 


(4.34) 

(4.35) 

(4.36) 

(4.37) 

(4.38) 

(4.39) 


We rewrite 


V h ’(p)q 


+ Pi9 


0 “Ms h-\\o + <*- l|V*?7_||o + ||V*r/+||o 


\jh' + {p+)q+ Q + Pig\\v+\\o- [p1p|t 2 I(p-) 2 
+ a _ ||V*? 7 _||q - [p] g \\g- - (t7-)IIo + \J h '-{P-)q- 


(4.40) 


+ (J+ IIV ||q , 
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where we use “(Pi)” to denote the average of rj± over T 2 . We first estimate the second line of 
(14.401) . We deduce from (14.381) . Cauchy’s inequality, and Lemma 13.101 that for any k > 0, 


[p] 5 |T 2 |(r?_) 2 = 


Ma 


|T 2 |(p+) 2 


di|T 2 |( r /+)+ / Q++ / $+ 

J q_|_ J r2_j_ 


Mff(/h) 2 |T 21 

(P + ) 2 


< 


1 


i + - fa+r + 

K 


2 , Mai i + «) 


Q+ 


|T 2 |(p+) 2 Un* 

+ Cy/£%mm{£Z,V"}. (4.41) 


By Holder’s inequality, we have 
2 

'J \l„ |2 




q+ I < 


h + (p+)\q+ [ 




i 


yjh' + (p+)q + 


Jn + h' + (p+) 

d 3 p + _ |T 2 |(p+ - pi) 


>n+ 9 

Now we choose the value of k via 

(1 + «)(P + -Pi) 


9 


y / h' + {p+)q+ 


(4.42) 


,+ 


= 1 <(=)> K = 


Pi 


Then we obtain 


[pi 5(1 + k) 

|T 2 |(p+) 2 


P 


9+ I < 


P^ - Pi 


(4.43) 


(1 + «)(P + -Pi) 

(p+) 2 

y/ti + {p+)q+ 


y / /i' + (p + )<7+ 


(4.44) 


and by Holder’s inequality, 

Mg(pim 




(p+) 2 

We deduce from (|4.41l) . (14.441) and (|4.45l) that 

2 

> ii - ,m 


1 3/91 |t 2 |, ^^Mapi 

—r—(p+) <—— 


(4.45) 


yJ h +(P+)q+ Q + Pi9h+\\l - M3|T 2 |(?7_) 2 

( sjh' + {p+)q+ q +pip||p+||o^ -Cy/S°min{eZ,V*} ( 4 ' 46 ) 
£ ll?+llo + IIP+llo - 

Next we use the estimate (I4.27P to obtain a bound for the first two terms on the third line 
of OH : 

l|V*p_||o - [pi 9 II 7 ?- ~ (p-)llo > (^- ~ <*c) ||V*t/_||o . (4.47) 

From (14.391) . Holder’s inequality, and Lemma 13.101 to obtain 

II 7 ?—Ho = II 7 ?- - fa-) llo + I t2 |(3-) 2 £ l|V*P-Ho + q-^j + 


< 


y/h'_(p-)q- + ||V*p_||o + y/£°Tam.{£°,T%). 


(4.48) 


Then we may combine (14.471) and (14.481) to estimate the full third line of (14.401) : 


v- llv*p_||o - [p! 9 II 7 ?- - ( 7 ?-)llo + \J h '-(p-)q- 


+ (7+ ||V*77+| 


> a + ||V*P+|| 2 + min{l,cr_ - u c } ||??_|| 2 - \/££min{££, V a n } (4.49) 
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We may then conclude from (14.401) . (14.461) . and (14.491) that 


V h '(p)<l Q + Pi9 ll*7+llo + a + ||V*?7+||o - [pj 9 11*7-115 + o ||V*r/_||o 
> ||g||o + min{l, a + } ||» 7 +||? + nun{l,<r_ - a c } ||?7-|| 2 - y/£° min{££,X>£}. 


(4.50) 


This is (14.33)) when a = 0. 

To derive (14.331) for a\ = a 2 = 0, ao > 0 we first take the time derivatives in (14.381) (14.391) . 
Then we may argue as above to derive the desired estimate. □ 


5. The evolution of energies controlling <9 3 q 

In this section we identify a dissipative structure for d 3 q and derive some energy-dissipation 
estimates. We again assume throughout this section that the solutions obey the estimate 
G% N (T) < 5, where 5 G (0,1) is given in Lemma fA.31 


5.1. Identifying the dissipative structure. Note that the energy evolution results presented 
in Sections I4.1H4.2I are not enough for us to get full energy estimates by applying Stokes regu¬ 
larity estimates as in the incompressible case [321 (33| . The problem is that we do not yet have 
control of diva. To control divu we appeal to a structure first exploited by Matsumura and 
Nishida [18J in the case in which p is a positive constant. We consider the following quantity, 
which is the material derivative of q in our coordinates: 

Q ■= d t q - Kd t 9d 3 q + UjAjkd k q = d t q - G 1 ’ 1 = - di v(pu) + G 1 ’ 2 . (5.1) 

From (13.241) we find that Q obeys the equations 

d 3 Q + p<9 3 (div u) = <9 3 G 1,2 - di v(d 3 pu) - d 3 pd 3 u 3 , 
pd t u 3 + pd 3 (h’(p)q) - pAu 3 - (p/3 + p')d 3 (divu) = Gj. 


We then eliminate d 33 u 3 from the equations (15.21) to obtain 


4^/3 + p! 

P)P 2 

- d t u 3 - 


(h'm) + Mh'm = 

dp/3 + p! 


dp/3 + p' 


d 3 G^ 2 + -G\ + d 3 h'(p)Q 


dp/3 + p ', 
h'(p)p 2 


(div(d 3 pu) + d 3 pd 3 u 3 ) + ^(d u u 3 + d 22 u 3 - d 31 u\ - d 32 u 2 ). 

P 


(5.3) 


P 2 


By the definition (15.11) of Q, we may view (15.31) as an evolution equation for d 3 q. This equation 
resembles the ODE dtf + f = g, up to some errors, and this ODE displays natural decay 
structure. We will extract this kind of structure for the more complicated equation (15.31) . 


5.2. Estimates. We now estimate the energy evolution of d 3 q at the 2N level. 

Proposition 5.1. Fix 0 < j < 2N — 1 and 0 < k < 4 N — 2j — 1. Then there exist universal 
constants A k ',j for 0 < k' < k so that 


dt 2 Xk ’’ j 

k'<k 


1+ 


k'<k 


< 

r^j 


d{ +1 u 

2 

+ 


4N—2j — l 


+ ^2 \^i N ~ 2j ~ k '~ l d 3 +l &i ( h'(p)q ) + | v m-2j-k'-i d k'+idiQ 

° k'<k 

v ^ N -2j-k-idiQ 2 | v* m - 2j - k 'di* 

° k'<k 

+ K(a, lp\)j£^V° 2N + K(a, M)^2 N^n + 2^n + K(a, [pIM+zKn. (5-4) 


2 

fc'+l 
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Proof. We first fix 0 < j < 2 N — 1 and then take 0 < k < 4 N — 2j — 1 and 0 < k' < k. Let 
a € N 2 so that |or| < 4 N — 2j — 1 — k!. Applying d a d 3 df to (15.31) and multiplying the resulting 
by d a d 3 +1 d{(h'(p)q) + d a d 3 +1 d{ ( h'{p)Q ) and then integrating over fi, we obtain 

I+ 11 +III = IV, 


where 


III = 


(h'(p)Q)) 

( tw' * (VfflS) , 

[ d°e$ ,+1 di(h'{p)q) 

Jn 


II = 


+ / d a d$' +1 di(ti(p)q)d a dl' +1 di {ti(p)Q) 

Jn 


and 


IV = £ {5“^' +1 ^(h'(p)g) (h'(p)Q)} 


x d-dk'Qj { 4 > l / : \ + I 1 ' d:iG V + + 


3 * 1 p 2 

4/x/3 + pi 


h'(p)p 


z^-d 3 h (p)Q - d t u 3 


(div(<9 3 pu) + d 3 pd 3 u 3 ) + ^(dnit 3 + «9 2 2^3 - d 3 i«i - <9 32 u 2 ) }> . (5.9) 

P 2 P J 

We will now estimate 1,11,111,1V. First, using the Cauchy-Schwarz inequality, we may 
easily estimate 


(5.5) 

(5.6) 

(5.7) 

(5.8) 


IV < 


d a d k 3 ' +1 di(h\p)q ) + £ I d a d k " +1 d{Q 


k"<k' 


dfG 1,2 


4N-2j 


+ 




4N-2J-1 


+ J] | d a dl"diQ 


k"<k' 


+ 


dj + \ 


4N-2j- 




k"<k' 


d a dl''d{u 


+ 


d a dl V*V3/u 


For the last term in III we recall the definition of Q from (15.11) in order to rewrite 

[ d a d k 3 ,+1 di(h'(p)q)d a d k 3 ' +1 di ( h\p)Q ) 

Jn 

= [ d a d k 3 +1 (i(h\p)q)d a d k 3 +1 4 (ti{p){d t q - G 1 ’ 1 )) 

Jn 

= \j t -J a ira^+ 1 a>(h'(p)q)d a ^+ 1 a{ (a'O)g 1 - 1 ). 

For II we estimate by expanding with Leibniz: 


II > 


4/r/3 + p' 


Jn h'(p)p 
For I, we have 

I > 


7.2 


k"<k' 


(5.10) 


(5.11) 


8 a d k ' +1 di (ti(p)Q) -C d a d k ' +1 d} (ti(p)Q) Y I d a dfd}Q 


o 

(5.12) 


J a rdf* 1 # ( ti{p)Q) d a d k 3 +i d{(h>(p) q ) 

d a d k 3 +1 di ( h\p)q ) Y, I d ° d l%Q 


-c 


k"<k‘ 


0 
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1 d f 4/i/3 + p' 


/ 


2 dt J n h'(p)p 2 

4^/3 + p' aa akp 


d a d k 3 ' +1 dt(h\p ) q ) 


j Q ^'(pjpi d a d3 +1 dKh'(p)q)d a d k 3 ' +1 di (h!(p)G 1,1 ) 


-C 


d a d k 3 +1 di(h'(p)q) J2 | d a d k "d{Q (5.13) 


k"<k’ 


Combining the estimates (15.101) (15.131) with (15.51) and applying Cauchy’s inequality in order 


to absorb the term 


d a d k ' +1 di ( h’(p)q) 


onto the left, we arrive at the inequality 


d_ 

dt 


'1 + 


4p/3 + p' k ' 


h’{p)p 2 


d a d k ' +1 diih\p)q) 

1 


+ 


d a d k ' +1 d{ ( h\p)q) 


+ 2 


d a d k 3 +i di(h\p)Q) < y 


k"<k’ 


d a d k &> t Q 


+ 


+ 


d{G 1 ’ 2 


i + d a d k ' +1 di(ti{p)q)d a d k ' +1 d{ ( h\p)G u ) 


AN—2j 


+ 


h'(p)p 2 

%G 2 2 


AN-2j-l 


+ 


d{ +1 u 


AN—2j—l 


+ 


V, 4N - 2j - k 'dju 


k'+1 


. (5.14) 


Owing to the Leibniz rule and the properties of p, we may estimate 


d k ' +1 d a d{Q " < ti(p)d k ' +1 d a diQ ' < d 3 +1 d a d J t (h'(p)Q) + Y 


k"<k' 


dfd a d}Q 


(5.15) 


Combining this with (I5.14|) and summing over all a with |a| < 41V — 2j — 1 — k', we deduce 
that there exist universal constants Pk\j,lk\j > 0 so that 


d_ 

dt 


'l + f +l di(h'(p)q) 


Pk' 


^AN^j-k'-lgk'+lgj (^) g ) 


+ Pk' 


^AN-V-k'-lgk’+lgjQ < 7fc/j ^2 \vi N - 2j - k '- 1 df9 J t Q 


k"<k' 


+ ,j Yj 

\a\<AN—2j—l—k' 




die 1 * 


AN-2j 


+ 


1 + 


&IG 2 


4/r/3 + p' 
h'(p)p 2 
2 

4V-2j-l 


0“9 3 fc,+1 ^(/i'(p)g)0“5 3 fc ' +1 ^(h / (p)G 1 ’ 1 ) 


+ 




4V-2j-l 


+ 


v ^4A-2 j-k'Qj u 


k'+l / 

(5.16) 


for every 0 < k' < k. We may then use Lemma IA.9I to deduce that there exist constants 
^k',j,a > 0 (depending on Pk',j,a and 7 k’,j,a as i n (1A.55|) ) and 


k,j ^k',j{lk',j + 1) 


(5.17) 


k'<k 


so that 


dt 2 

k'<k 


1 + 


+ S (h'(p)g) + II v AN-2j- k '-i d k'+idiQ 


k'<k 


k’<k 
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\a\<4N—2j—1—k' 
k'<k 


in 


1 + 


4^/3 + p'\ „ Q „ fc ' 


h'(p)p 2 


d a d k 3 +1 &>{h\p)q)d a d$ +1 &{ {h!{p)G 1,1 ) 


+ Afcj 


( 

dje 1 ’ 2 

2 

+ 

&>G 2 

2 

+ 

d{ +1 u 

V 


AN—2j 


AN—2j—l 



+ Afcj 


4IV-2J-1 

v iN-2j-k-lQjQ 2 + ^2 \W* m ~ 2j ~ k 'd{ 


k'<k 


2 

fc'+l 


(5.18) 


Finally, we will estimate the nonlinear terms in the right hand side of (15.181) . We use the 
estimates (13.361) of Lemma 13.31 to estimate, for 0 < j < 2N — 1, 




+ 


4N—2j 

Then we use Lemma 13.71 to bound 


%G 2 


< 


41V—2j—1 


-^(<T [pl)^2JV^2Ar + ^N+2^N- (5.19) 


Ki E / (i + a“# +1 #(ft'G9«)a"# +1 a? (fe'(p)G 1 ' 1 ) 

|a|<4JV-2j-l-fc' ^ ' 

k'<k 


< 


K(<r, lpi)V^N\/ £ 2N V 2N + £n+2^2N (5-20) 


Plugging the nonlinear estimates (15.191) and (15.201) into (15.181) then yields (15.41) . 


□ 


We then record a similar result at the IV + 2 level. 

Proposition 5.2. Fix 0 < j < N + 1 and 0 < k < 2(N + 2) — 2j — 1. TTien t/iere exist universal 
constants \k',j f or 0 < k' < k so that 


d_ 

dt 


fc'<fc 



4/r/3 + // 2(AT+2)-2j-l-fc' 

/!'(p)p 2 V * 


% +1 di(h'(p)q) 


2 

0 


+ E ll v 


2(JV+2)-2j-l-fe' 


aj' +1 8f (a'(p)«) 


fc'<fc 




fc'<fc 


< 

re 



2 

+ 

4N—2j — l 


v 2(N+2)-2j-l-k^ Q 


2 

0 


+ x; \\vi (N+2) - 2j - k 'diu 

k'<k 


2 

fc'+l 


+ K(a, Ip})V£2N V n+2- 


(5.21) 


Proof. The proof proceeds along the same lines as that of Proposition 15. II except that we instead 
use the estimates (13.411) of Lemma 13.41 □ 


Remark 5.3. Propositions \5.1\ and \5.2\ 'provide two important bits of control: energy estimates 
of d 3 q and dissipation estimates for d^Q. These are crucial for improving the horizontal energy 
and dissipation estimates derived in the previous section into the full ones in later sections, 
respectively. 


6. Combined energy evolution estimates 

Now we combine our previous estimates with the elliptic regularity theory of a particular 
Stokes problem in order to derive an intermediate energy-dissipation estimate. We again assume 
throughout this section that the solutions obey the estimate Gf N (T) < 5, where 6 £ (0,1) is 
given in Lemma lA.31 
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6.1. The Stokes problem. We first derive elliptic estimates. We deduce from (13.2411 that 
di v(pu) = G 1,2 - Q, 


— 4An — ^ 4~ - V div u + V (h'(p)q) = —G 2 — dtu. 
P P P 


( 6 . 1 ) 


Direct calculations give the form of the Stokes problem we shall use: 
-pA + V (h'(p)q) = ^zG 2 - d t u - p ^2<9 3 d 3 u + d 33 Q 

^/ 3 ±i4 V (^(G^-Q-d^us) 


+ : 


div ( - J = — (G 1,2 - Q - 2d 3 pu 3 ) 


u = u 


in D± 

in f2± 
on d£l±. 


( 6 . 2 ) 


We now prove the Stokes estimates at the 2N level. 

Lemma 6.1. Fix 0 < j < 2N — 1. Then for any 1 < k < 4 N — 2j, 
2 


v 4V-2j-fcm 


u 


+ 




fc+i 

< 

r^i 

dl +l u 

2 

+ 

AN—2j — 1 


W* 4N ~ 2j ~ k di ( h\p)q) 


2 

fc-1 


V: 


AN—2j—k 


dfQ 


+ ||^* Ar ' u ||i + K(a, IpDF^n^n + ^n+2-^2N- 


(6.3) 


Proof. We first fix 0 < j < 2N — 1 and then take 1 < k < AN — 2j. Let a £ ft 2 so that 
|a| < AN — 2j — k. Applying d a d{ to the equations (16.211 in respectively, and then applying 
the elliptic estimates of Lemma IA.11I with r = k' + 1 > 2 for any 1 < k' < k, by the trace 
theory, we obtain 


d a d{u 


fc'+l 


+ 


< 


d a d{G 2 


Vd a d{ (h'(p)q) 

+ 


2 

< 

d a dt( 4) 

2 

+ 

k'—l 

\pj 

fc'+l 


Nd a d{ ( h\p)q ) 


k'-l 


k'-l 


d a di +1 


u 


+ 


d a dfG 1,2 


+ 


d°d]Q 


k’-l 


+ 


+ 


d a d 3 t u 


k’ 


d a d 3 u 


< 


d a d 3 G 1 ’ 2 


+ 


d a d{G 2 


k -1 


+ 


d a d? +1 u 


H k ’+ 1/2 (E) 
2 


fc-1 


+ 


d a dj Q 


+ 


d a d}u 


(6.4) 


k' 


+ 


V/d a dju 


H i/2(E) 


< 


&> t G 1 ’ 2 

2 

+ 

9{g 2 

2 

+ 

d{ +1 u 

2 

+ 


AN-2j 


4AT—2j—1 


4IV—2j —1 


V: 


AN-2 j-kgjQ 


+ ||V^u||J + 


d a d 3 t u 


A simple induction based on the above yields that 


d a cP t u 


fc+i 


+ 


Vfrtfq 


+ 


t { 

d} +1 u 


fc-1 

2 


&t G 1 ’ 2 


AN-2j 


+ 


dlG 2 


+ 


N, 4N - 2j - k djQ 


m-2j-i 

Finally, we use the estimates (I3.36P of Lemma 13.31 to have 


4IV—2j —1 

+ ||V^|li 


(6.5) 


die 1 ’ 2 


AN—2j 


+ 


%G 2 


< 


AN—2j—l 


K(a, 


2N'E > 2N + £n+2^2N 


We then sum over such |a| < AN — 2j — k to conclude (16.31) . 
We then record a similar result at the N + 2 level. 


( 6 . 6 ) 

□ 
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Lemma 6.2. Fix 0 < j < N + 1. Then for any 1 < k < 2 (N + 2) — 2j, 


v ^2(IV+2)-2j-fe^ u 


2(IV+2)+l 


+ 




k -1 


< 

r^j 


dj +1 i 


2(N+2)-2j-l 


+ 


V^N+2)-2j-kgjQ 


+ 


V 


2(AT+2) 


+ K(a, \p\)£-2n'D < n+2- 


(6.7) 


Proof. The proof proceeds similarly as Proposition 16.11 by using instead the estimates (13.411) of 
Lemma l3~H □ 

6.2. Synthesis. We will now chain the energy evolution estimates of Section [4] with the d 3 q 
estimate of Section [5] and the estimates of Lemmas 16.1116.21 The full dissipation estimates of u 
will be obtained, and also some estimates of q will be improved along the way. To do so, we 
first introduce some notation. For n = 2 N or n = N + 2 we write 


C= / p\Vl n u\ 2 +ti(p)\Vl n q\ 2 + I p 1 dVfr /+ r + a + |V*Vf7 ?+ 


rz2ri 1 2 


w2 ra„ 1 2 


in 


+ J ~ [p! 9 |V* n p_| 2 + cr_ |V*V* n p_j 2 

+ [ (p(J-l)K + q + d 3 pe)J\d?u\ 2 + h'(p)(J-l)\d?q\ 2 (6.8) 
Jn 


and 


= 


[ ^ |B 0 V+x| 2 + p'|divV 2n - 12 

Jn 2 


u\ 


(6.9) 


for the various terms appearing in Propositions I4.1l and 14.31 Similarly, for n = 2N or n = IV+ 2 
and integers 0 < j < n — 1 and 0 < k < 2n — 2j — 1 we write 

2 


:= J] Afc/j 
k ’=0 


i/ 1 + 


( 6 . 10 ) 


with the constants \k',j the same as in Proposition 15.II in the case n = 2N and as in Proposition 
5.21 in the case n = N + 2. We also write 


* j n ■= E {h'(p)q) + £ IIVV 


k '=0 


7 2n-2 


k'=l 


fc'-l 


+ llvr*c|i;+ e , ( 6.ii) 


k '=0 


where Q is as defined in (15.11) . In addition, we introduce the following intermediate energies: 


£n := HvNIIo + EIK? 
1=0 


2n—2j 


+ 77(— [pi) 


E*> 2ri _2 ■ + rninjl. cj} E|| 
1=0 n 3 3=0 


2n—2j+l 


+ 'W([p])min{l,(7 + ,a_ - ^c} E ^ 

l=o 


2ra-2j+l 


V n : = E |^t l 
1=0 


n— 1 


271—2^+1 


+ E v ^+(p)g) 


l=o 


2n—2j—1 


( 6 . 12 ) 


(6.13) 


The rest of the section is devoted to the derivation of the energy bounds for and V n based on 
the energy evolutions for (£((,£>((,21+,23+. First we present the main energy evolution result 
at the 2 N level with the improved dissipation energy T> 2 N- 
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Proposition 6.3. Let ® 2 /v be defined as above and let T> 2 N be as defined by (16.131) . 

There exist universal constants 7 2 jv > 0 and fi 2 N;j,k > 0 for j = 0,... , 2N — 1, k = 0,..., AN — 
2j — 1 so that 


d_ 

dt 



2N-1 4N—2j—l 


+ 5Z Pwj,k Zl 2N 

j =0 k =0 


+ T2n'E > 2N 


< K(a, [p]) y/SfaVS N + K(a, M)^V" N £° N+2 T 2N + K{a, [p])^ +2 ^ 2 n, (6-14) 


Proof. First, we sum the result of Proposition 14.11 with j = 2N, with the result of Proposition 
EM for all a € N 1+2 with ao < 2N — 1 and |a| < 41V; this yields the estimate 

J t <%N + ®2N < K{a, M)V^ V 2 N + K(T m)^N^N + 2^N, (6.15) 

where <£ 2N an d ® 2 n are as defined in (16.81) and (16.91) . Note that, owing to the Korn inequality 
of Proposition IA.81 we may bound 


|V^«||J <2> 


2N- (6-16) 

Next, we recall the notation Q in (15.11) . We may use the estimates (|3.36D of Lemma [3731 to 
obtain the bound 


V! 


r Q||o < II Vf div(^) 11 0 +11 V^G 1 - 2 11 o < \\vi N u\\l+K(aM)^ a 2 N +£ a N + 2 ^ 2 N. (6.17) 
Then for 0 < j < 2N — 1 and 0 < k < AN — 2j — 1, we may combine the results of Proposition 
5.11 and Lemma 16.11 (summed over 1 < k! < k) with (16.171) to see that 


d 


dt 2N w *o 0M ^ 


k 

2N rsj 


dj +1 u 


4N—2j—l 

+ K(a , Ip}) 


+ y hy- 2j ~ k, cp f Q 


k '=1 


k' 


+ ||vr«||J 


£f N V°2N + \/7?9V^’v-i_9-7 r 2Ar £n4-2^ 


2N c N+2- r 2N' ~r c N+2- r 2N 


, (6.18) 


where we have written 21^ and as ai A6.101) and (16.111) and employed Lemma f6.II to control 


the term 


S/^N-2j-k'Qj u 


k '+1 


in the right hand side of (15.41) . Notice that if we write 


y k N ■= iivrsiio+ y^-^diQ 


2 

fc+1 


then we particularly have 


y k N < i|vrs|io+E 


k '=0 


— -°2 N' 


(6.19) 


( 6 . 20 ) 


Let 0 < j < 2N — 1 and 0 < k < AN — 2j — 1. We sum the estimates (|6.151) and (|6.181) : 
employing (16.201) in the resulting estimate, we deduce that there exists a universal constant 
C\ > 0 such that 


n) + ^>2N + Cif>2N ~ H^ 7 * 4 ^ 2j k 


k '=o 


k' 


+ 


d! +1 i 


4N—2j—l 


+ * u||i + A (a, [pj) k/ £f N V 2N + J'£ > 2 n£n+2'£‘ 2N + £-n+2^2N (6-21) 


for all such j, k. 

Now, for fixed 0 < j < 2N — 1, owing to the definition (16.191) of P\ 2 n > we ma y y i ew A6.21B as 
a sequence of estimates indexed by 0 < k < AN — 2j — 1, of the form considered in Lemma lA.91 
Applying the lemma, we deduce that there exist universal constants C2-jk > 0 so that 


4iV—2j —1 

J t E 

k =0 


4N—2j—l 

+ %2n) + 

k =0 


2) 


2N 


+ CiS) j ’ k 


2N 


< 


;+ 


2 

4N—2j—l 
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+ 


< 


V/ N ~ 2j diQ 


d] +1 u 


+ IE u ||x +(o') M) V ^2N^2N + v'^ > 2N^N+2^2N + ^N+2-^2N 


1 + ®2iV + K{ a i [pD V^N^N + ^ a 2N^N+2^N + E% +2 F 2N j (6.22) 


4N—2j — 

where the second estimate follows from (16.1711 and (16.161) . 
It is easy to see, using the definitions of 'S^N an d ^ 2 N■ 


v^||? + ||vrQ||o+ d}Q 


2 

4N—2j 


4N—2j—l 

< E (^2N + ClSi 


k= 0 


j,k 

2N 


(6.23) 


Adding 


di +1 


u 4Ar _ 2 ._ 1 + K{(T, M)£Zn V 2N + ^n+ 2 ^ 2 N to both sides, and using Lemma 16.1 
(with k = AN — 2j ), we then have that 


d{u 


+ V& t {h\p) q ) 

AN—2j + l 
AN—2j—l 

~ E (®2N + C^2n) + 

k =0 


4Af-2j-l 


di 


+ 1 , 


u 


AN—2j—A 


+ K(a, [p] )^2N^2N + ^n+2-^2N- (6.24) 


Hence there is a universal constant C^-j > 0 so that 
4AT—2j—1 


"j7 E C 2-j,k{^2N + ^2n) + C ‘A'J 


dt 


k =0 


dfu 


AN-2j+l 


+ 


Vdi(h'(p)q) 


4N—2j—l 


< 

r^j 


dj +1 u 


47V _ 2j ._i + ®2iv +-^( a ! M) ^n + ^ 2N S^ 2N + S^ 2 n (6.25) 
for all 0 < j < 2N — 1. 

Counting backward from 2N — 1 to 0, we may view (16.2511 as a sequence of inequalities of 
the form considered in Lemma IA.91 Applying the lemma, we find that there exist universal 
constants C^ ; j > 0 so that 


d 


2N-1 AN—2J—1 


2N-1 


E E ( &2N + 2 I 2 A) + E ^3'j 

j =0 k =0 j =0 

< lln2AL.ll 2 


d{u 


AN-2j+l 


+ 


V%{ti{p)q) 


4N—2j— 1 


u||x + [pD \f^N^2N A- \JT>2 n^n+2 -^2 N + £n+2 ^ 2 N 

~ ® 2 AT + -^(‘L [Pi) a/ ^2N^2N + '^ > 2 n^N+2-^2N + £n+2^2N j (6.26) 
where the last inequality follows from (16.161) . 

Let C 5 > 0 denote the universal constant appearing on the right side of the last inequality 
in (16.2611 . We multiply (16.151) by 2 C 5 and add the resulting inequality to (16.2611 . This results in 
the bound 


d_ 

dt 


2N—1 4N—2j—l 


+ E E 

C'i;jC 2 -j t k‘^2N 

j =0 k =0 




2N-1 / 


+ C^ 2 N + 

E C 3; 3 ( 

d{u 


3=0 V 



AN-2j+l 


+ 


Vdi(h'{p)q) 


4N—2j—\ 


< 

r^j 


K{a, [pj) + J^n + 2^n + £ a N+2 F 2N , (6.27) 


where 


2IV—1 AN-2j -1 

Cq = 2 C 5 + E E r ':/' 2 ,M- 

j =0 k =0 


(6.28) 
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Finally, we estimate 


2 N 


2N—1 


E* 

3=0 


U 


4N-2j+l 


+ E \ vd i( h \p)<i) 

3=0 

2N—1 

<C 5 V° 2N + c 3;j 

3=0 


4N—2j—\ 


dt u 


4N-2j+l 


+ 


Vdi(h’(p)q) 


4N—2j—l 


(6.29) 


This and (16.271) imply (16.141) upon dividing by Cq and renaming the universal constants. □ 

Next we record a similar result at the N + 2 level. 

Proposition 6.4. There exist universal constants 'Jn +2 > 0 and £isr+ 2 ;j,k > 0 for j = 0,..., N+ 
1, k = 0,..., 2 (N + 2) — 2j — 1 so that 


d_ 

dt 


7V+1 2(AT+2)-2j-l \ 

&N+2 + E E. PN+2-,j,k^N + 2 J + lN+2^ ) N+2 % K(<J, [p])\/ ^f^N+2 • (6.30) 

j =0 fc =0 J 

Proof. The proof proceeds as Proposition 16.31 using instead Propositions 14.2114.4115.21 Lemma 
16.21 and the estimates (13.41 [) of Lemma 13.41 □ 

Next, we show that, up to some error terms, (££ is comparable to £° for both n = 2 N and 
n = N + 2. 

Lemma 6.5. Let £f be defined by (16.121) . be defined by (16.81) . and 6 e defined by (16.101) . 
Further let ^ 2 N;j,k an d PN+ 2 -,j,k be the constants appearing in Propositions E3 and \6-4\ respec¬ 
tively, for appropriate integers j,k. Then 


,k 
*2IV 


2JV—1 m-2j-l 

£%N — yj£%N min {^2IV, ^2 at} ~ ^2IV + E E ^2N-j,k^21 

j= 0 /c=0 

^ -^(<L [pJ)£flV + vElV min {^2IV) d^2N } (6.31) 

and 

7V+1 2(N+2)—2j — l 

£n+2 ~ a/ £f N min{^+2- ffiv+ 2 } ~ ^A+2 + E E PN+2;j,k^N + 2 

J=0 fc=0 

< iL(fj, [/0])<?Ar_|_2 + y^ m in{^ +2 ,^ +2 }. (6.32) 

Proof. We will only prove (16.311) : the proof of (16.321) follows from a similar argument. Let us 
compactly write 

27V-1 4IV—2j —1 

Z:= E E (6-33) 

j= 0 /c=0 

First note that 


2IV-14IV-2J-1 2N-14N-2j-l k 

Z =E E ~ E E E||vf'-«-‘ , - i ^' +i a?(ft'(p) S ) 

j =0 / c =0 j =0 / c =0 / c '=0 

27V—1 47V—2^—1 27V—1 

= E E |vf , -«-*- i %^(k'(p) g ) | x x) ||m'('>'(?)9) 

j=0 k= 0 " j=0 

Since 

1 


. (6.34) 

4IV—2j—1 


d 3 d}q = —- d 3 d 3 t {h'(p)q ) - d 3 (/i'(p))<9^J . (6.35) 

and h'(p) is smooth on [— 6 , 0 ] and [ 0 , £] and bounded below, we may estimate 

' 2 < 2 + 11^-^111 ( 6 - 36 ) 


d 3 d}q 

2 

< 

0 ~ 


2 

+ 

0 

dlq 

2 

<2( + 

0 ~ 

^q 
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and similarly 


d 3 d{q 

2 

< 

2 

d 3 &t (h'(p)q) 

2 

+ 

i 

d}q 

2 

< z + 

2 

Vt d}q 

2 

+ 

0 

d 3 djq 


2 

2—1 

<2 + 11^-^111 + 


d 3 d}q (6.37) 

2—1 


for i = 1,..., 4 N — 2j — 1. This constitutes a sequence of estimates of the form considered in 
Lemma IA.9I (setting Xi = 0 there to remove the time derivative terms). Applying the lemma, 
we find that 

AN—2j—l 


d 3 d 3 t q 


AN—2j— 


i < d3d t q 


i =0 


< Z + ||V 4 A -1 I 


On the other hand, 


dlq 


AN—2j 

so summing (16.3811 yields the estimate 

2N 

1=0 


<||V^||o + 


d 3 cP t q 


4 AT—21—1 


AN—2j 


<z+m N q \\l 


(6.38) 


(6.39) 


(6.40) 


Next we note that Proposition 14.61 implies that 

*Zn > IIvr?Ho + ||V^llo + n- [pD [IIV^pIIo + * IIV*V^||o 

+ ^(M) min{l, cr_|_, <r_ - a c } ||V^7/||J - y/£!famin{eg N , V% N } 

+ [ ( p(J - 1 )K + q + d 3 p0)J\d 2 t N u\ 2 + h'(p)(J - 1) \d 2N q\ 2 . (6.41) 
Jo. 

We may easily estimate 


[ (p(J - l ) K + q + d 3 p9)J\df N u\ 2 + h'(p)(J - 1 ) \ d 2N q^ 2 
Ja 


min {£% n , V% n }, (6.42) 


l v *Mlo+Hl v * v *Nlo ~ 1 p 7 ? 2n _ 2 - +min { 1,CJ }X] I H 7 ? 

3=0 n ° 3=0 


2n—2j+l 


and 


3=0 


2n-2j+l 


(6.43) 


(6.44) 


Plugging (I6.42|) - (|6.44|) into (]6.41 1) then yields the improved bound 


«2N>||V^||o + ||Vr«||o + ^(- [pD 


n 

2 || 

2 

3= 0 

+ mini Lex} > \\dip 

2n—2j 1 4—f || * 

1=0 

2n—2j+l 


+ 'W([pl)min{l,( 7 + ,c 7 _ -a c }J2\\dlrj 

3=0 


2n-2j+l 


~ V^n^o.^n^n}- (6-45) 


To conclude the proof of the first inequality in (16.3111 . we sum (16.4011 and (16.4511 and employ 
(16.331) . The second inequality in (16.3111 follows easily from (16.341) . (16.421) . and the definition of 
(££ in (16.811 and £% N in (16.121) . 

□ 
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7. Comparison estimates 

In this section, we shall show that £% * s comparable to £° and that T>° is comparable to T> n , 
for both n = 2 N and n = N + 2. This will be achieved by deriving further elliptic estimates. 

We again assume throughout this section that the solutions obey the estimate G% N (T) < 5 , 
where 5 € (0,1) is given in Lemma I A. 31 


7.1. Energy comparison. We begin with our comparison result for the energy. 

Theorem 7.1. Let £° and £° be as defined in (12.151) and (16.121) . respectively. It holds that 

&2N ~ [p]) (z%N + (£%.n) 2 + £n+2^2N + (£%.n) 3//2 ) (7-1) 


and 


£ 


N+2 £ K{a, [pj) 


(7.2) 


£n+ 2 + ^2N^N+2 + {&N+2) 

Proof. We let n denote either 2 N or TV + 2 throughout the proof, and we compactly write 

11/2' (7-3) 


w, = ||v 2 ” _ 2 g 1 ||5 + ||v 2 ” _ 2 g 2 || 2 + ||v“ _ 2 g 3 || 2 /2 + ||v 2 ” _ 1 g 1112 


Note that by the definitions of £f in (16.121) . we have 

+ hi{— [pi 


C > K 


u 


lo + Z)|^ 

7=0 


2 

2n—2j 


3=0 


2n—2j 


+ min 




7=o 


2n—27+1 


+ W([pJ) 


min{l, cr_|_, <r_ — a c } Ell din 
7=0 


2n-2j-\-l 


- (O s/2 


(7.4) 


We first estimate dju for j = 0,... ,n — 1. The key is to use the elliptic regularity theory of 
the following two-phase Lame system derived from (|3.24l) : 


—fiAu — (/i/3 + p')V divu = G 2 — pdtu — pV (h'(p)q) in fl 

—§(u)e 3 = ( -P'(p)q + p\gr)+ - <r+A* 77 _|_)e 3 + G^ on E + 

- [§(«)] e 3 = {lP'(p)qj + [pi gi 7 - + cr„ A*r/_)e 3 - G 3 _ on E_ 

[tij =0 on 

U- = 0 on E/,. 


(7.5) 


We let j = 0,..., n — 1 and then apply c£ t to the problem (17.51) and use the elliptic estimates of 
Lemma lA. 101 with r = 2n — 2j > 2, by (17.31) (17.41) and trace theory to obtain 


d%u 


< 


2n—2j 


d{G 2 


+ 


2n—2j—2 

dlv" 


+ 


di +1 


u 


2n—27—3/2 


~\~ O 


2n — 27 — 2 
2 


+ 


d{q 


2n—2j—l 


+ 


dig 


H 2n-2j-3/2(£) 


< 


d{ +1 u 


dj.V 

2 


2n—27+1/2 


+ 


diG 3 


2n—2/—3/2 


+ iL(a,[pl)fe+W n + (^) 3/2 )- (7-6) 

2n-2(j+l) V / 

Using a simple induction based on the estimate (17.61) . utilizing the ||<9 ”u||q estimate in (17.41) for 
the base case, we easily deduce that for j = 0 ,..., n, 

2 


dju 


< 


2n—2j 


K(a, m)[^ + yVn + (S‘ 


(7.7) 


We then estimate djq and d^p for j = 1, • •., n to get an improvement. By the first equation 
of (I3.24p . using the estimates (17.71) and (17.31) . we have that for j = 1,..., n, 


dig 

2 

< 

d{ l u 

2 

+ 

di~ l G l 


2n—27+1 ~ 


2n—27+2 



2n—2j+l 
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d{ 1 u 

2 

+ 

^" 1 G 1 


2n— 2(j- 1) 



2n—2(j — 1)- 

Now by the kinematic boundary condition 

d t r] = u 3 + G 4 on X), 

we have that for j = 1 ,,n, by the trace theory, m and (EH), 


1 <^,[p])(^+Wn + (^) 3/2 )- (7.8) 

(7.9) 


dj/n 


2n—2j+3/2 


< 


' u 'i 


< 

r^j 


di S 


2n-2(j-l) 


+ 


^2™-2j+3/2( S ) 

i,2 

ST G 4 


+ 


^ _1 G 4 


2n-2(/-l)-l/2 


2n-2j+3/2 

<K(a, [pl)fe + W n + (^ 


Consequently, summing the estimates 117.41) , 117.71) , (17.81) and (17.101) , we conclude 

£n<K(a, [pJ)fe+Wn + (^ 


(7.10) 

(7.11) 

Setting n = 2iV in (17.111) . and using ([3.351) of Lemma [3731 to bound W 21 V < iv(cr, [p])[(^ 27 v) 2 + 
£n+ 2 -F‘ 2 n], we obtain (17.11) : setting n = N + 2 in (17.111) . and using (13.401) of Lemma [374] to bound 
wv+ 2 < K(a, iP\)£? 2 N £ N+ 2 : we obtain E2J). □ 

7.2. Dissipation comparison. Next we consider a similar result for the dissipation. To begin 
we prove a lemma involving norms of q. 


Lemma 7.2. Let T> n be as defined in (16.131) . Then 


12 n 


+ 


IS- 


Proof. We compute 


Vq= W) l V (h'(f>)q)-V(ti(p))q] 


Since 6/(p) is smooth on [—6,0] and [0,^] and bounded below, we may estimate 


ll^gllo < ||V(/l'(p)g)|| 2 + Iloilo < Pn + 


and similarly 


l|Vg|| 2 <||V(/ l / (p) (? )|| 2 + || (? || 2 <Pn + 


lo + HVglllr 


(7.12) 

(7.13) 

(7.14) 

(7.15) 


for i = 1,..., 2rt — 1. This constitutes a sequence of estimates of the form considered in Lemma 
IA.9I (setting X,; = 0 there to remove the time derivative terms). Applying the lemma, we find 
that 

2n—1 

l|Vg|| 2 n _i< E ]|Vg||?<£n + 

i=0 


IS. 


(7.16) 


which yields the desired estimate upon adding ||g|| 0 


□ 


to both sides. 

Now we can record the dissipation comparison results. 

Theorem 7.3. Let V £ and V n be as defined in (12.161) and ([6.131) . respectively. It holds that 

(7.17) 


^2N ~ Tv (<7, [p]) ( T>2N + £2n'E > 2N + + £n+2^2N 


and 


V 


N + 2<K(a, Ip})(Vn+2 


+ £2n'£ > N+2 + \! £ 


'N+2 U N+2 


(7.18) 


Proof. We again let n denote either 2 N or N + 2 throughout the proof, and we compactly write 


T n = ||V 2ri - 1 G 1 || 2 + llV 


IIV^G 3 " 2 


+iiv^- i gi; /2 


+ \\V 2n - 2 d t G 


1/2 

4 " 2 +ct 2 ||V 2n G 4112 


(7.19) 


11/2 1 II * II1/2 • 

We now estimate the remaining parts of T> n not contained in T>f, that is, to estimate rj and 
improve the estimates of q. We divide the proof into several steps. 

Step 1 — q estimates 
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We first notice that by the first equation of (|3.24j) . 


Wlbn-l — 11 ^ 112n + ||^ ILn-l ~ 


and for 2 < j < n + 1 , 
dtQ 


2 

< 

d{ l u 

2 

+ 

dl^G 1 

2n-2j+2 

2n—2j+3 


< 

di _1 u 

2 

+ 

di^G 1 


2n—2j+3 



2n-2j+2 
2 

2n-2j+2 


(7.20) 


(7.21) 


+ Tn- 


We may sum these estimates with the estimate of Lemma 17.21 to see that 

n+2 

+ Tra + 


lL + ll^llL-l + \\^tQ 

3 =2 


2n-2j+2 


2 

0 • 


(7.22) 


Step 2 - dir] estimates 

We now derive estimates for time derivatives of r/. For the term djp for j > 2 we use the 
kinematic boundary condition 

dtf] = 1 x 3 + G 4 on E. (7.23) 

Indeed, for j = 2,..., n + 1, by the trace theory and (17 .19 [) . we have 

2 


d} v 

2 

< 

d{ 1 u 3 

2 

+ 

^“ 1 G 4 


2n-2j+5/2 ~ 


#2n-2j+5/2( E ) 



2n-2j+5/2 


< 


d{ 1 u 

2 

+ 



2n—2(j—1)+1 



2n-2(j-l)+l/2 

For the term dtp, we again use (17.231) . the trace theory and (17.191) to find 


+ Tn- 


(7.24) 


V 2 II^IIL+1/2 + ll^hllL-1/2 ~ (! + (j2 ) Il“3||^2n+1/2( E ) + cr 2 ||G 4 || 2n+ i/2 + IIII2 ti—1/2 

< (1 + max{cdj_, cr 2 _ }) \\u\\ 2 2n+1 + y n < K{a } \fi\)V n + y n . (7.25) 

Hence 


a 11 ^7112n.+l/2 + \\dt r l\\2n-l/2 


n+1 


+ ^2\\ d tv 

3 =2 


2 

2n-2j+b/2 


<K(a, M)v n + y n . 


(7.26) 


Step 3 - V*r/ estimates 

For the remaining p term without temporal derivatives we use the dynamic boundary condi¬ 
tion 


- cr + A*p + + pigr] + = P[(pi)q + - 2 p + d 3 u 3t+ - p! + divu + - G| + on E + (7.27) 

and 

- cr_A*?/_ - [p]] = - [P'(p)g] + 2 {pd 3 u 3 J + [p'divtt] - on E_. (7.28) 

Notice that at this point we do not have any bound of q on the boundary E, but we have 
bounded V(h'(p)q) in 9. As such, we first apply V* to (17.271) and (17.281) . We use Lemma fA. 121 
on <1727)1 : when [p] < 0 we also use Lemma [A. 121 on (|7.28D . but when [p] > 0 we instead use 
(I4.24p of Lemma 14.51 This, trace theory, and (17.191) then provide the estimate 


■H(M) min{cr 2 , (cr_ - a c f} ||'V*77 || 2n+1/2 + U{- Jp]) (V ||V*p || 2n+1/2 


+ ||V* 


12n—3/2 


< 

r^j 


V*G 


3||2 


I V*(/l (p)<Z) ||#2n-3/2( S ) + ||V*Vn|| H 2n-3/2( E ) + || v *^3 11 2n-3/2 


< 


||V(h / (p)g )||^_ 1 + ||u || 2 n+1 + \\G 3 \\ 2 2 n _ 1 / 2 <V n + y n , (7.29) 


where P = X(0,oo) denotes the Heaviside function. 


Step 4 - Recovering 


and 
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Next we seek to control ||(?||q and || 77 1|q. To this end, we recall the conservations of mass 
(I4.38p - (|4.39l) and the definition of <f> in (I3.65|) . Standard arguments allow us to hnd w = w±, a 
solution to the following problems: 

divu;_)_ = q + + 4> + in f2 + , u; 3i+ = — pip_|_ on E + ,ws t+ = —p + p_ on S_ (7.30) 

and 

divum = g_ + <f>_ in W 3 - = —p~p- on w- = 0 on £*,. (7-31) 

Note that the equalities (|4.38l) and (|4.39l) provide the necessary compatibility conditions for 
the solvability of the problems (17.301) and (|7.31D . respectively. Moreover, we have the following 
estimates 

IMIi ^ lkll 0 + Iloilo + IMi/ 2 - ( 7 - 32 ) 

Now we take the dot product of the second equation of (|3.24l) with ( p)~ 1 w , we obtain 

f X7(h'(p)q) ■ w — div§(u) • (p)~ 1 w = f ((p ) _1 G 2 — dtu)w . (7.33) 

J Q J Q 

Using the equations for w = w±, we have 

[ V(/i' + (p+)< 7 +) ■ w+ - divS + (u+) • (p + )~ 1 w + = [ {p + )~ l (P' + {p + )q+I - S+(u+))e 3 • w + 

J f2-|_ J E+ 

- [ {P+)~ 1 (K(P+) < 1 + 1 ~^+( u +)) e z-w+- [ h' + {p + )q + div w + - S + {u + ) :X7({p + y 1 w + ) 
J e_ J 

= ~ f V+ [{P[(p+)q+ I - §+(n+))e 3 ] • e 3 - (pi ) -1 Y {§ + (u + )) i 3 w i!+ 

• 7s + 2=1,2 

+ [ V- [( P' + {p+)q + I - § + (u+))e 3 ] • e 3 - (p +) _1 Y (S + (u + )) a w it+ 

• 7s - 2=1,2 

- [ h' + {p+)q + {q + + <f>+) - §+(«+) : V((p + ) _ 1 u; + ), 

J r2_|_ 

(7.34) 

and similarly, 

/ V (h'_(p-)q-) ■ w _ - div§_(u_) • (p_) _ 1 u;_ = - [ p- [( P'_{p-)q-I - §_(u_))e 3 ] • e 3 
Jft- Je_ 

~Y [ (P-) _ 1 (S-(“-))i 3 Wi,- - [ ti_(p_)q_(q_ + $_) -S_(u_) : V((p_) _ 1 u;_). 

2=1,2 • 7s - 


(7.35) 


Collecting ()7.33f) (17.35D . and the boundary conditions of (I3.24D . we get 

[ {{p)~ l G 2 — d t u)w = — [ ti(p)q(q + <f>)- S(u) : V((p)~V) 

J £2 J Q 

- [ V+ [(P'+(p+)q+I - S+(w+))e3] • e 3 ~ Y (Pi) _1 ( S +( u +))*3^,h 

• 7s + 2 = 1,2 

+ f V- [T’ , (p)g/-S(w))e 3 ] -e 3 - [(p) _1 ( s N)*3^i] 

2 = 1,2 

= ~ [ h / (p)g( (? + $)-§(u):V((p)- 1 u;) 

Jn 

- / ^+(PiOT+ - cr+A*r? + + G| + ) - V (pi) -1 (§>+(u + ))i 3 Wi :+ 

■ /s + 2=1,2 

+ f V-(IpI 9V- + °-^*P- - G|-) - Y [(^) _1 ( S (“))i3^i] ■ 

2 = 1,2 


(7.36) 
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We then further deduce from the above that 

I h'{p)\q\ 2 + f Pl g\r) + \ 2 - f {p} g\r)-\ 2 + f a |V*p| 2 
Jn Jt,+ J e_ Jt, 

= ~ [ {{pr l G 2 - dtu)w - [ ti(p)q<f> -S(u): V((p) _1 u;) - [ pG\ 

J q J Q J s 

- [ (/h) -1 ( s +( u +))w+ - [ [(p) _1 (s(«))i3^] ■ 

i= 1,2 "' s + t=l,2’' E - 


(7.37) 


Using Proposition 14.61 trace theory, and Cauchy’s inequality, we deduce from the above that 
for every e > 0 


I l + W-i-M) Ikllo + <r\\v*v\\l 


kilo + 'H(M) mi n{ 1 ,o-+,cr_ ~a c } 

~ Iloilo IMIo + Iloilo Iloilo + Ikllo ll $ llo + Nil 
+ Ikllo N 3 ||o + l|Vn||^ (E ) NIIhO(E) + 


1C 


1 


< 


C 


o + Nil + Iloilo + N 3 ||o + Iloilo 


+ £ 


13 + 


+ ni Pl)y/£ZK, (7-38) 


where again U = X(o, oo) ■ 

Step 5 - Handling 1)7.3811 with cases 

To proceed we must break to cases depending on the sign of [pj. If [p] > 0 then we employ 
(17.3211 and Lemma 13.101 to successively estimate 

£ (iMIi + Ikllo + IMIS) < CEdkllS + IMII? + IWlS) < CedkllS + 11,11? + s/SSK)- (7.39) 

We then choose e > 0 according to 

Ce = ^min{l,cr + , (cm - a c )}, (7.40) 

so that we may absorb ||(/||q + ||p|| 2 onto the left in (17.381) . This yields the estimate (again 
estimating 4> with Lemma 13.101) 


Ikllo + mm{ 1 ,a + ,a- - a c } ||? 7 || 2 

< K(a, [p]) (j|<kklo + ||u||| + II^Ho + ||G' 3 || 0 + y/S^V^j 

< K(a, [p]) (V n + y n + y[E°v£} , (7.41) 

where the constant K(cr, [p]) is as in (|2.40|) when [p] > 0. 

On the other hand, if [p] < 0, then we use (17.321) and Lemma f3. 101 fto estimate 4>) to bound 


w 


l? + 


I II2 , 

kilo + 


10 j ^ Ce 


112 . 
?llo + 


<F 


<Ce(\\q\\ 2 0 + 


l 2 + 

.lo + 

| 2 + 
lo + 


IV* 


13 + 

vfo 


v*k 
+ y/e*K)- 


1 11 112 

In this case we choose e > 0 so that Ce = which allows us to absorb ||g|| 0 + 
left of (17.3811 . From this we deduce (again using Lemma f3.101) that 


(7.42) 


onto the 


I 112 . 

kilo + 


+ <7 IIV* 


< 

r^j 


2 | II 112 | || v'Y 

n + u o + G 


'o + \\° 3 


where in the last estimate we have used (17.291) to estimate ||V*p|| 0 . 
Combining (17.41|) and (|7.43l) with (17. 2911 . we find that 


0 + llv*klo + VW 

<Vn+yn + V£ZVn, ( 7 - 43 ) 

2 


Ikllo + 77([p]) min{l, o + , cm - a c , <j|, (<r_ - cr c ) 2 } |kllL+ 3/2 
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+ W) (lMlL-1/2 + min{l,cr 2 } hllL+3/2) 

< K(a, [p]) (v n + y n + y/SZV*) . (7.44) 

Step 6 - Completion 

Now we sum (17.221) and (17.261) and add to the resulting inequality a positive multiple of 
(17.441) . where the multiplier is chosen large enough to absorb onto the left the term ||(/||q on the 
right side of the sum of (17.221) . This results in the estimate 


n+2 


MIL + \\ d tq\\L-i + ^2 |p q 

3 =2 


2 n—2j+2 


+'H(M) min{l, CT+, cj_ - (j c , crL (o'- -cr c ) 2 } IMlL+3/2 


+ %(~ M) 


-1/2 + min{l, a 2 } \\r]\\ 2 2n+3/2 j + ll^llL-1/2 + l|dMlL+i/2 

<K(a,lp})(v n +y n + ^£ZV"). (7.45) 


n+1 

+ 5Z|| 9 f r? 

3 =2 


2n-2j+5/2 


When n = 2 N we use (13.361) of Lemma HOI to estimate TW K(a, IpD^^n^n + £n+ 2 ^n- 
Then we obtain (17.171) from (17.451) by adding T> 2 N to both sides and recalling that by definition, 
(16.131) , T> 2 n controls all of the u terms in • When n = IV + 2 we use (13.411) of Lemma 13.41 
to estimate 3V+2 K (<r, Then we obtain (17.181) from (|7.45l) by adding T>n +2 to 

both sides. 

□ 


8. A PRIORI ESTIMATES 

We are now ready to derive the global-in-time bounds and decay of high order energy £% N 
and £^+2 based on the previous estimates on various energies. 

Recall that Qf N is defined by (12.241) . 


8.1. Boundedness at the 2 N level. In this subsection, we shall show the boundedness at 
the 2 N level. We first recall J~ 2 N defined in (12.171) . We will make use of the following transport 
estimates for p 2 N- 


Proposition 8.1. There exists a universal constant 0 < ho < 1 so that if (t) < 6 < 5q, then 

T 2 jv(r) < J2At(0) + r f Vf N (8.1) 

Jo 


for 0 <r <t. 

Proof. The estimate is recorded in Proposition 7.2 of [8J. The proof is based on the transport 
theory of the kinematic boundary condition. □ 


We now present the a priori energy bounds for £f N , P^ad an d J~ 2 N- 

Proposition 8.2. There exists a constant I\(cr, |[p]) of the form (12.401) and a universal constant 
5o > 0 so that if 0 < 5 < 5o/K(a, [p]) and Qf N {t) < 5, then 

sup £f N {r)+ I Vf N {r)dr+ sup < K(a, [p] )£% N (0) + F 2 n (0)■ (8.2) 

0 <r<t Jo 0<r<t (T + H 

Proof. Throughout this proof let us denote by £f N the quantity appearing in the time derivative 
in (16.141) of Proposition 16.31 By assumption we have that sup 0<r<t £f N (r) < 5. In what follows 
we will restrict the size of 5 in order to prove (18.21) . We assume initially that <5 < 1 so that 
(£% N ) Pl < {^ 2 n) P2 T Pi > P 2 - This allows us to simplify many of the subsequent estimates by 
retaining only the lowest power of £f N appearing in inequalities. We also assume that 6 is as 
small as in Lemma lA.31 
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We may integrate the estimates (|6.14l) of Proposition 16.31 to see that 


£ 


2 N 


(*) + /^<^( 0 ) + ™ j sj SZn'DZn + \J' 


' £> 2N £ N+2-^ 2N + £ N+2^2N 


N+2 J 


(8.3) 


According to the first inequality of (16.311) in Lemma 16.51 combined with estimate (17.11) of 
Theorem o we have that 


£ 2N < K ( a i [pD £ 2N + (^2v) 3 ^ 2 + ^N+2^2N 

Estimate (17.171) of Theorem 17.31 yields the bound 


^2N ~ IpD ^2 N + V^2N^2N + £n+2-^2N 


(8.4) 


(8.5) 


Additionally, the second inequality of (16.311) in Lemma 16.51 together with the trivial bound 
£- 2 n — £ 2 N> implies that 

£ 2N ^ K (t [p])£ 2 V (8-6) 

We may assume, by taking the maximum, that the constants K(a, [p]) > 1 appearing in (18.31) 
(18.61) are identical. 

Let us assume that V5K(a, [p]]) is sufficiently small to absorb the second terms on the right 
sides of (18.41) and ([8.51) onto the left with a factor of 1/2. Then (18.41) and (18.51) improve to 


£ 2N < -K (<r, [/?]) £%n + £n+2^2N 


and 


^2JV ~ (O', [/?]) \P2N + £ N+2^2n\ 

We then plug (|8.6I) - (|8.8D into (|8.3I) to see that 

£ 2Nit) + [ £> 2N [Pl)^27v(0) + K (°5 lPi) £ N+2(t)J 7 2N{t) 

JO 


(8.7) 

( 8 . 8 ) 


T K(a, M) 2 J \/ £ 2 n' £> 2N + ^om£m^2N + £/V-i_?-£"2. 


'2N u 2N 4- \j i^2N c N+2- r 2N ~r Z N +2- r 2 N 

For any e > 0 we may apply Cauchy’s inequality to bound 


K(a,M) 2 I* 

Jo 


nr c-n t- ^ , K MPi ) 4 /*« 

U 2N t 'N+2 J ~ 2 N — £ j U 2N H ^ / t -N+2- r 2N- 


10 


^2 


(8.9) 


( 8 . 10 ) 


io 


Taking e sufficiently small, we may absorb the /(j V% N term onto the left in (18.91) , resulting 
the bound 

£ 2n(^) + [ T>2N AT(fJ, [p])^2A(0) + lpi) £ N+2(t)J 7 2N(t) 


m 


£ 


+ K(a,M) 2 [ \^£^V^ N + {l + K{aM?)£N + 2^2N\. (8.11) 

Jo 1 

The decay of £^+ 2 ^) guaranteed by the bound on combined with the bound of 

Proposition 18.11 easily imply that 

(t)F 2 N(t) < <5(1 + t)~ 4N+8 T 2 N (t ) < <LF 2A r(0) + <5 f* T> 2 N {r)dr 

Jo 

<VS^ 2 n( 0 ) + VS f T> 2 N {r)dr (8.12) 
Jo 

t £ a N+ 2 {r)T 2 N{r)dr < 5T 2N ( 0) + 5 t V a 2 N (r)dr < VSF 2N { 0) + V5 f V% N {r)dr. (8.13) 
Jo Jo Jo 


and 
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We plug (|8.12|) and (|8.13l) into (18. lip and then further assume that Vd(K(a, [p]]) + K(<j, [p]) 2 + 
K(a, [p]) 4 ) is sufficiently small to absorb all of the resulting JqU^n t erms on the right. This 
results in the estimate 

&2N r(t) + f ^wi^dr < K(a, [p])£f N (0) + ^NiO). (8-14) 

Jo 

Then the estimate (18.21) follows from (|8.14j) and (18.11) . □ 

8.2. Decay at the N + 2 level. We next show the decay at the N + 2 level. The first step is 
an interpolation argument that allows us to control £%t, 2 in terms of , 2 and £f N - 

Lemma 8.3. Let £%j, 2 and £f N be as defined in (12.151) and 'D'fi +2 ^ e as defined in (12.161) . and 
let 

0 = St! 6(0,1). 


If [p]| > 0, a + > 0, and cr_ > a c , then 


'N+2 


< max < 1 


AN — 7 


min{l, ct|, (T — a c } 


min{l,fj + ,cr_ - a c ,cr+, (<r_ - a c ) 2 } 


V 


N+2- 


If M < 0 an d o+,o- > 0, then 


£n+2 — max { 1) —: 


minjl, cr_|_, <r_} 


V c ,! 


If M <0, then 
In either case, 


min{l, 0 - 2 ,0-2 }J ^+2' 
£n+2 ~ (I^N +2 ) ^ i.^2N ) 1 9 ■ 


£ n +2 ~ M)(^ > Af+2) 6l (^2Ar) 


1-6 


(8.15) 

(8.16) 

(8.17) 

(8.18) 
(8.19) 


Proof. The estimates (18.161) and ()8. 171) follow directly from the definitions. 
Suppose then that [pj < 0. First note that we may trivially estimate 


N+2 

J2 \H u 


1=0 


N+2 


2(N+2)-2j 


+ lkll2(AI+2) + E 


N+2 

< \\di u 

1=0 


1=1 


2{N+2)-2j+l 


N+2 


2(N+2)—2j+l 


+ 


1=1 


2(Y+2)—21+3/2 


+ IMl2(V+2) + ll^t"" 2 


2(V+2)—1 


< m+2ns™) 1 - 6 ( 8 . 20 ) 


2(N+2)-2j+2 


N +3 

+ II 2(Y+2)—1/2 + S \\^ q \ 

1=2 

since all the terms appearing in the middle can be controlled by both , 2 and £f N ■ It remains 
only to estimate 

ll»lll2(JV+2) + mini 1 , o'} \\v\\1(n+2)+i ■ ( 8 -21) 

To handle these remaining terms we must use Sobolev interpolation. Indeed, we first have that 

ll»lll(K+2) < IMl“tf+2)-l/2 IIXIIIT"’ < 

where 6 is as in (18.151) . To handle the second term, we interpolate twice: 
min{l, a} || 27 II 2 ( 7 V+ 2 )+i < mm{l,u} IN| 2(A r +2 ) +3 / 2 IMl 2 (JV+ 2 )+i /2 


( 8 . 22 ) 


< min{l, a} |M| 2(JV+2)+3/2 IMI 2( v+ 2) _i /2 HUn ' 

= (min{l,a 2 }||r ? || 2 (7V+2)+3/2 ) ' ( M\\n+ 2 )-i /2 

< N+2) e (£m) l ~ e , 


2 \ 4iV-7 

4 N , 


(8.23) 
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where the last inequality holds because 

1 AN — 9 

2 + 81V- 14 


4IV — 8 
AN - 7' 


(8.24) 


Then (|8. 18fl follows from (18.201) . (18.221) . and (18.2311 . 

Now when [p] < 0, (18.191) follows from (18.181) since K(a, [p]) = 1 is a constant of the form 
(12.401) . When [pj > 0 we chain the trivial bound T>f^ +2 < + 2 )°{^ 2 n) 19 with (18.161) and 

note that 


K[a, IpJ) 


max 


min{ 1, , a- — a c } 

' min{l, cr + , cr_ - a c , <j\, (<r_ 



(8.25) 


is a constant of the form (12.401) . 

□ 


Next we deduce algebraic decay of £fj +2 - 

Proposition 8.4. There exists a constant K(a, [p]) of the form, (12.401) and a universal constant 
do > 0 so that if 0 < 6 < 60 /K(a, [p]) and Gf N (t) < 5, then 

sup (1 + r) iN ~ & £ a N+ 2 {r) < K(<t, [pl)^(0) + F 2 N (0). (8.26) 

0<r<t 


Proof. First we note that again sup 0<r<t £f N {r) < 5, and that we will restrict 6 to prove the 
desired result. Let us initially assume that 5 is as small as in Proposition ^. 21 which in particular 
means that S < 1. In this proof we will also need to explicitly keep track of various universal 
constants Cj > 0. Throughout the proof we will write £^, 2 for the term appearing with the 
time derivative in (16.301) of Proposition 16.41 and we will write 


1 

4N-8' 


(8.27) 


We begin by enumerating various estimates proved previously. From Proposition 16.41 we have 
that 

—£n + 2 + Ci^n+2 ^ C 2 K(a, lpj)^/£f N T > N + 2- (8.28) 

By assuming that VSK(a , [fpll) is sufficiently small, we may deduce from (17.21) of Theorem ITT 
and ([USD of Theorem O that 


£n+2 < C 3 K(aM)£N+2 (8-29) 

and 

V a N+2 <C 4 K(a,mm +2 . (8.30) 

Further restricting 5 in this manner, we may use the first estimate in (16.321) of Lemma 16.51 and 
(18.291) to see that 

0 < ~£n + 2 < £n +2 _ V£^n^n+2 < C 5 &N +2 (8.31) 

Similarly, the second estimate in (16.321) of Lemma 16.51 and the trivial estimate £pj +2 < £n +2 
imply that 

£n +2 <C 6 K(<t,IpI)£n +2 - ( 8 - 32 ) 

Next, from (18.191) of Lemma 18.31 we know that 

£n+ 2 < C 7 K(a, [pi) (T>n +2 ) ^ {£ 2N ) ■ (8.33) 

Finally, since 5 is as small as in Proposition 18.21 we know that 


sup £ 2N (r ) < C% (. K[a , [pJ)£ 2 V(0) + ^2Jv(0)) := M 0 . (8.34) 

0 <r<t 

As before, we are free to assume that the constants K(a, [p]) > 1 appearing in ()8.28D — (18.30D 
and (|8.32I) - (I8.34I) are all the same. We may also assume, increasing the stated constants if need 
be, that 


C 4 > C\ and Cq,C 7 > 1. 


(8.35) 
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Again restricting 5, we may guarantee that C 2 C 4 K ( a , [[/o ]]) 2 \/£ 2N < Ci/2, which allows us 
to chain together (|8.28D and (18.301) to find that 


—£^+2 + -^' D ‘ n +2 - 0 ) 


which in turn implies, by way of (18.301) . that 

Lfr + 

1 .C-N+2 + 


dt 


2C 4 iv (c, 


:'E > N+2 < 0. 


(8.36) 


(8.37) 


Chaining together (18.3211 — (18.3411 and using the positivity of £y +2 guaranteed by (18.311) shows 
that 

(C 6 C 7 K(a, lpj) 2 ) 1+s M s 0 ( ^+ 2)1+S - P ^+ 2 ’ (8 ' 38) 

where 0 < s < 1 is as in (18.271) . We may then combine (18.371) and (18.381) to see that 

d_ 

dt 

where we have written 

Ci 


—Sn+2 + Z(£n+ 2 ) S ^ 


(8.39) 


2 1 2C4(C 6 Cry+*K(a, m) 3+2s M* 0 ' (8 ' 40) 

Now we view (18.391) as a differential inequality for £^ +2 , which is positive by virtue of (18.311) . 
We may integrate (18.391) to find that for any 0 < r < t, 

&N+ 2(0) 


^v_L9(r) < — _ 

+ - [1 + 5 Z(£S +2 (0))*r]l/. 

From (18.321) and (18.341) we know that 

&N+2 (0) < Cs£" n+2 { 0) < C 6 K(a, [p])f 2 ^(0) < C 6 K(a, lpj)M 0 , 
which implies that 

sCi{£" n+2 ( 0 )) s sCi(C 6 K(a,lp})M 0 y 


(8.41) 

(8.42) 


sZ(£" n+2 (0)) s = 


< 


2C 4 (C 6 C7) 1+s Ar((j, [p]) 3+2s A4o - 2C A {C,C 7 y^K{aMf +2s M s Q 

= 2C A C & {C 7 y+°K{aMY +s - 1 ’ (8 ' 43) 

Here in the last inequality we have used (|8.35|) and the fact that K(a, [pj) > 1 and s < 1. A 
simple computation reveals that for 0 < M < 1 

(1 + r) 1//s 1 


r>o (1 + Mr) l / S M l / S ' 


(8.44) 


From this and (18.411) (|8.43|) we then know that 
(l + r)V^ +2 (r)< (! + 0 1/s 


[1 + sZ(£° N+2 (0)yr]Vs-' N+2 '-' - s s Z s £" n+2 { 0) 

(2Ci) 1 / s (C 6 C 7 ) 1 / s + 1 K(a, 1pI) 3 / s + 2 

{sCif/s 

This immediately yields (18.261) upon recalling the definition of Mo hi (18.341) . 


( 0 )< 


'^+ 2 ( 0 ) 


Mo. (8.45) 


□ 


Finally, we deduce exponential decay of £y , 2 in the case with surface tension. 

Proposition 8.5. Suppose that either [p] < 0 and <t + ,ct_ > 0 or else [pj > 0 and a + > 0, 
ct_ > <t c . Define M(cr, [[p]) > 0 according to 

min{l, <r+, <r — cr c } 


M(a, | 


= max < 1, 


min{l,(T + ,c7_ - a c ,al, (cr_ - a c ) 2 } 


(8.46) 
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if [pj > 0 and 


M(a,\p\) = max 


min{l, o + , ov} } 
’ min{l. J 


(8.47) 


if M < 0. 

There exists a constant K(a, [p]) of the form (|2.40[) and a universal constant 5o > 0 so that 
if 0 < 6 < Sq/K( a, [p]) and G 2N (t) < 5, then 


sup exp 

o <r<t 


( K{a , \pj)M(a, [p]) 


^N+2( r ) ~ K(a, [p])£jv-+2(0). 


(8.48) 


Proof. Assume that 5 is as small as in Propositions 18.21 and 18.41 Arguing as in Proposition 18.41 
(and renaming the constants), we know that 


+ k£m v ’™ - 0 

and 

0 < jj-jj^V+2 ^ &N +2 ^ C 3 K(a, |[pl)^ + 2- 

Also, from (18.161) and (18.171) of Lemma 18.31 we know that 

&N +2 — c a m(o, mm +2 , 

where M(a, [p]) is as dehned above. 

Combining (18.491) — (18.511) leads to the differential inequality 


d_ 

dt 


£ 


N+2 


+ 


Cl 


C 3 C 4 K(a, [pJ) 2 M(cr, [pj 


-£ 


N +2 4 ; 


< 0. 


Upon integrating and again using (18.501) . we then find that for 0 < r < t, 

rC\ 


£n +2 ( r ) 


C 2 K(a, [ 

This immediately yields (I8.48p . 


^ &N+2 (r) < C 3 K(a, [p])^ +2 (0) exp 


C 3 C A K(a, [pfl) 2 M(u, 


(8.49) 

(8.50) 

(8.51) 

(8.52) 

(8.53) 

□ 


8.3. A priori estimates for G 2N . Now we combine the results of Propositions 18. ll 18.21 and 
18.41 into a single a priori estimate involving the functional Qf N , as defined in (|2.24l) . 

Theorem 8.6. There exists a constant K(cr, [p]) of the form (12.401) and a universal constant 
So > 0 so that if 0 < 5 < 5o/K(cr, [pj) and Gf N (t) < S, then 

G 2N {t) < K(a, M)£% n (0) + T 2N (0). (8.54) 

Proof. Let <5o > 0 be as small as in Propositions 18.1118.21 and 18.41 Then the estimates of each 
Proposition hold, and they may be summed to deduce (18.541) . □ 


9. Proof of Main results 

9.1. Proof of Theorem l2.31 Before presenting the proof of Theorem l2.31 we record a technical 
result that we will use to transition from the bounds provided by the local well-posedness result, 
Theorem 12.21 to bounds on Q 2N , as defined by (12.241) . 

Proposition 9.1. Suppose that N > 3. There exists a universal constant C > 0 with the 
following two properties. 

First, if0<T, then we have the estimate 

G 2 jvC^ 1 ) — SU P ^ 2 n(^) + [ + SU P J~ 2 N{f) + C(1 + T) iN 8 sup £f (9-1) 

0 < t<T Jo 0 < t<T 0 < t<T 

Second, if 0 < T\ < T 2 then we have the estimate 
G 2 n(T 2 ) < CQf N {T\ ) + sup £ 2 m(t) + f T> 2N (t)dt 

Ti<t<T 2 .7Ti 
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+ 7T—TFT SU P J 7 2 N(t) + C(T 2 -T 1 )\l + T 2 ) 4N - 8 sup SS N (t). (9.2) 

(I + JlJ T!<t<T 2 T\<t<T 2 

Proof. Although the exact form of the energy and dissipation terms is slightly different, the 
result follows from the same argument used in the proof of Proposition 9.1 of [8]. As such, we 
omit further detail. □ 

With this result in hand we may now turn to one of our two main results. 

Proof of Theorem \2.S[ The structure of the estimates in our local well-posedness Theorem 12.21 
our a priori estimates Theorem 18.61 and our Proposition 19.11 are essentially the same as those 
found in corresponding results from the study of the incompressible viscous surface wave problem 
carried out in [8]. As such, we may argue as in the proof of Theorem 1.3 of [8] to prove (I2.28p . 
Rather than replicate the argument here, we will provide only a sketch; we refer to [8] for a 
more detailed version of the argument. 

Step 1 - Local theory 

Set 5 > 0 to be the smaller of the So constants appearing in Theorem 18.61 and Proposition 18.51 
divided by the larger of the K(a , [pj) constants appearing there. By choosing k < So, where 
<5 0 comes from Theorem 12.21 we may guarantee that Theorem 12.21 provides a unique solution 
on an interval [0,T] with T < 1, satisfying the estimates (j2.181) . We then use estimate (19.111 of 
Proposition 19.11 and further restrict k to guarantee that Qf N (T) < 5. By the choice of 5, the 
estimates of Theorem 18.61 and Proposition 18.51 hold. 

Step 2 - Global existence 

We define 

T*(k) = sup{T > 0 | for every choice of initial data satisfying the compatibility 

conditions and + F2jv(0) < n there exists a unique solution on [0, T] 

that achieves the data and satisfies Qf N C T)<6 }. (9.3) 

By Step 1, we know that T*(/t) > 0, provided n is small enough. Clearly, T* is non-increasing, 
and so it suffices to establish that there is a kq > 0 so that T*(«o) = +oo. 

To prove this we argue by contradiction, showing that T*(ko) < +oo leads to a contradiction 
if no > 0 is chosen appropriately. More precisely, we choose a 0 < T\ < T*(ko) and then use 
Theorem 18.61 to guarantee that the hypotheses of the local existence theorems are verified by 
(q(Ti),u(Ti),r](Ti)). We then apply the local theory to construct a solution on [Ti,T 2 ], By 
using estimate (EH) of Proposition 19.11 and choosing T\ sufficiently close to T*(ko), we may 
guarantee that T*(«o) < T 2 and Qf N (T 2 ) < 5. This contradicts the definition of T*(ko)- So, 
T*(ko) = Too. It’s clear from this analysis that l/«o is of the form K(a, [p]). 

Step 3 - Proof of (12.281) 

To deduce the bound ([2.281) we simply apply Theorem 18.61 on the interval [0, oo). This is 
possible since the previous step guarantees that Qf N (+°o) < 5. 

Step 4 - Proof of (12.291) 

Since 5 was chosen as small as in Proposition 18.51 the hypotheses of this Proposition are 
satisfied in the cases ([2.2611 and (12.271) . As such, the estimate (18.481) holds, from which (12.291) 
easily follows. 

□ 


9.2. Proof of Theorem 12.91 We can now establish the vanishing surface tension limit. 

Proof of Theorem \2.9\ For the sake of brevity we will only present a sketch of the proof. 

The convergence of the data occurs at sufficiently high regularity to pass to the limit in the 
compatibility conditions, which shows that (<7o> uo,rjo) satisfy these with a± = 0. Then the data 
convergence and the bound (12.281) show that 

lim sup G 2 n{°°) ^ &2N + -^2Ar(0). 

(<T_|_ ,<T— )—^0 


(9.4) 
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These bounds, when combined with the usual embedding and interpolation results, show that 
the solutions (q a , u a , rj a ) converge strongly to a triple ( q,u,r 7 ) in a regularity class high enough 
to pass to the limit in (12.101) . The limiting ( q,u,r 7 ) then solve (12.101) with a± = 0, and hence 
agree with the solution obtained from Theorem 12.31 in this case. □ 


Appendix A. Analytic tools 


A.l. Poisson extensions. We will now define the appropriate Poisson integrals that allow us 
to extend r)±, defined on the surfaces £±, to functions defined on fi, with “good” boundedness. 

Suppose that E + = T 2 x {£}, where T 2 := ( 2ttLiT ) x ( 211 L 2 T ). We define the Poisson integral 
in T 2 x (— oo,l) by 


where for £ G (L 1 1 Z) x (L 2 1 Z) we have written 


m 



fix') 


p-ii-x' 

- -- dx'. 

2it\/ L1L2 


(A.2) 


Here ” stands for extending downward and stands for extending at x 3 = £, etc. It is 
well-known that V- (_ : H S (T , + ) —>• fp+V 2 (T 2 x (— 00 , £)) is a bounded linear operator for s > 0 . 
However, if restricted to the domain fl, we can have the following improvements. 


Lemma A.l. Let V- if be the Poisson integral of a function f that is either in or 

Pf q_ 1 / 2 (S + ) for q € N = {0,1,2,...}, where we have written H S (T , + ) for the homogeneous 
Sobolev space of order s. Then 

W q V-,if\\ 0 < ||/||^- 1/2(t2) and ||V?P _,,/|| 0 < \\f\\% q{T2) . (A.3) 

Proof. See Lemma A.3 of [ 8 ]. □ 


We extend r/ + to be defined on H by 

f}+(x',x 3 ) = V + ig + (x',x 3 ) := V-^rj+ix', x 3 ), for x 3 < £. (A.4) 

Then Lemma I A. 1 1 implies in particular that if rj + G fif s_ 1 // 2 (S + ) for s > 0, then f) + G H s {Pi). 
Similarly, for £_ = T 2 x {0} we define the Poisson integral in T 2 x (— 00 , 0) by 

(A'5) 


It is clear that has the same regularity properties as V-j. This allows us to extend t/_ to 
be defined on H_. However, we do not extend ??_ to the upper domain by the reflection 
since this will result in the discontinuity of the partial derivatives in x 3 of the extension. For 
our purposes, we instead to do the extension through the following. Let 0 < Ao < Ai < • • • < 
A m < 00 for m G N and define the (m + 1) x (m + 1) Vandermonde matrix V(Ao, Ai,..., A m ) 
by V(Ao, Ai,..., Xm)ij = (—A j) 1 for i,j = 0,...,m. It is well-known that the Vandermonde 
matrices are invertible, so we are free to let a = (ao, or, ..., a m ) T be the solution to 

V(Ao, Ai, • • • , A m) & — Qm , Qm — (1, 1, . . . , 1) . (A. 6 ) 

Now we define the specialized Poisson integral in T 2 x (0, oo) by 


V+flfix) 


E 


e j P x ' 

27T\/ L 1 L 2 


j =0 


(A.7) 


It is easy to check that, due to (IA. 6 I) . d l 3 V +} of(x', 0) = d l 3 V- t of(x', 0) for all 0 < l < m and 
hence 

d a V+ t of(x', 0) = d a V- : of(x', 0), Va G N 3 with 0 < |a| < m. (A. 8 ) 
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These facts allow us to extend r/_ to be defined on Ll by 


r]-(x',x 3 ) = V-r)_(x',x 3 ) 


'P+, oV-(x',x 3 ), x 3 >0 
P-, oV-(x\x 3 ), x 3 < 0 . 


(A.9) 


It is clear now that if Tj- G for 0 < s < m, then G H s ( fi). Since we will only 

work with s lying in a finite interval, we may assume that m is sufficiently large in (1A.6D for 
r/_ G H s (n) for all s in the interval. 


A.2. Estimates of Sobolev norms. We will need some estimates of the product of functions 
in Sobolev spaces. 


Lemma A.2. Let U denote a domain either of the form or of the form E-j-. 

(1) Let 0 < r < si < S 2 be such that si > n/ 2. Let f G H Sl (U), g G H S2 (U). Then 
fg G H r {U) and 

ll/S'll/D- ~ ll/ll/pu • (A.10) 

(2) Let 0 < r < si < S 2 be such that S 2 > r + n/2. Let f G H Sl {U), g G H S2 (U). Then 
fg G H r (U) and 

ll/sll/p- ~ II/IIiY s i llslltf-a • (A. 11 ) 

Proof. These results are standard and may be derived, for example, by use of the Fourier 
characterization of the H s spaces and extensions. □ 


A.3. Coefficient estimates. Here we are concerned with how the size of rj can control the 
“geometric” terms that appear in the equations. 

Lemma A. 3. There exists a universal 0 < 5 < 1 so that if IMI 5/2 < <5, then 


\J l|lL°°(n) + 


ll-AA-1 
II K 


l°°(o) + ll^llL oo (n) — 2’ 

1 


lL°°(r) '^11-^- l|lz,°°(r) — 2 ’ an d 


(A.12) 




+ 


L°°(n) 


< 1 . 


Also, the map 0 defined by ( 12 .4p is a diffeomorphism. 

Proof. The estimate (IA.12D is guaranteed by Lemma 2.4 of [BJ. 

A.4. Korn inequality. Consider the following operators acting on functions u : M n D U 


□ 


rp n 2 div u 

D u = Vu + S/u 1 and D°u = Du- 1, 


n 


or in components 


D mj = diUj + djUi and B%i 7 - = diUj + djUi -^ 6{j. 


n 


Note that tr(D°u) = 0. As such, the operator D° is referred to as the “deviatoric part of the 
symmetric gradient.” 

In what follows we will compute the kernel of D and D° in H 1 ( Q). In doing so we will actually 
compute the kernel on C 3 (f 2 ), but the results hold as well in via an approximation 

argument. 


Lemma A.4. 

ker(D) = {u(x) = a + Ax \ a G M n , A = —A T } (A.13) 

Proof. We have that Du = 0 if and only if diUj = —djUi for all i,j. Then 

didkUj — dkdjUi — Oj { fl/U/.) — c)jcf u/,. — 3, ( v,j ) — QjDigij , 
and hence V 2 u = 0. Then u(x) = a + Ax for some a G M n and A G M nxn . But then 

0 = Du = A + A t => A = -A t . (A.14) 

So, if u G ker(D) then u(x) = a + Ax for A antisymmetric. The converse also clearly holds. □ 
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~ C | u is holomorphic}. 


(A.15) 


We may also compute the kernel of B°. 

Lemma A.5. If n = 2 then 

ker(B°) = {u : C ~ R 2 D li — i 

If n > 3 then 

ker(B°) = {rt(x) = a + Ax + 7 x + (b ■ x)x — b\x\ 2 /2 | 7 € R, a, b € R n , A = — A T }. (A.16) 

Proof. Consider first the case n = 2. Then B°rt = 0 if and only if 

2d\ui diU 2 + d 2 u{\ _ fd\ui + d 2 u 2 0 


2d 2 u 2 




0 


d\u 2 + d 2 u 2 , 


K diu 2 + d 2 ui 
which holds if and only if 

d\U\ = d 2 u 2 and d 2 u\ = —d\u 2 , 

which are the Cauchy-Riemann equations in Q. Regarding R 2 ~ C then 
holomorphic in 12, which is ()A.15D . 

Now consider the case n > 3. Throughout the rest of the proof we do 

summation convention. Note first that if B°u = 0 then 

„ „ 2 div u 1 . 

2diUi = Brie* • e* =-re* • e, = — divrt 

n n 


for each i = 1 ,,n. This means that 


1 


d\U\ = d 2 u 2 = • • • = d n u n = — div u := M. 

n 

This allows us to rewrite B°ti = 0 in components as 

diUj = —djUi + 2 MSij. 

Next we compute 

did k Uj = d k (-djUi + 2M5ij) = -djd k Ui + 2 d k M5ij 
= —dj(—diU k + 2 M5ik) + 2d k M8 t j 
= didju k + 2 d k MSij - 2djM8 ik 
= di(-d k Uj + 2M5j k ) + 2d k M5ij - 2djMS ik 
= -did k Uj + 2d k MSij - 2djM8 ik + 2 diM5 jk , 

and hence 

did k Uj = d k M5ij - djMS ik + diM5 jk . 

Applying dt to (I A. 231) . we find that 

did k deuj = d k diM5ij - djdiM5 ik + did £ M8 jk . 

On the other hand, we can switch the index from k to I in (IA.23D to see that 

didkUj = diM5ij — djMbu + diM5jg, 
which reveals, upon applying d k , that 


(A.17) 

(A.18) 
shows that u is 

not employ the 
(A.19) 

(A.20) 
(A.21) 


(A.22) 


(A.23) 
(A.24) 
(A.25) 


did k deuj = d k diM8ij - djd k M5 ie + d t d k M A 7 . (A.26) 

Now we equate (IA.24I) and (1 A.26I) to see that 

d k diM5ij - djd(_M5i k + did(,M5 jk = d k diM8ij - djd k MSi£ + d % d k M <5^, (A.27) 

which allows us to cancel the d k dnM8ij terms to deduce that 

did e M5 jk - djdiMS ik = did k M5ji - djd k M5i £ . (A.28) 

Multiplying by 5j k and summing over j, k then reveals that 

(n - 2 )did e M = -AMfe, (A.29) 

from which we deduce, by multiplying by 5i£ and summing over i,£, that 

(n - 2) AM = -nAM. 


(A.30) 
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Hence AM = 0, and upon replacing in (1A.29D and using the fact that n > 3, we find that 

V 2 M = 0 =>• M(x) = 7 + b ■ x for some 7 €E R n , b € M n . (A.31) 


Returning to (IA.24D . we also find that 

V 3 u = 0. 

From (IA.31D and ()A.20D we deduce that 


(A.32) 


x 2 

Ui(x) = gi{xi) + 7 Xi + (b ■ Xi)xi + , (A.33) 

where Xi = x — x&i, and gi is (by virtue of (jA.321) 1 a quadratic polynomial in x t . i.e. in all of 
the variables except x*. 

Since whenever i / j we have djUi = — diUj , we may use (IA.33D to compute 


djgi{xi ) + Xj 6 j = - digj(xj ) - 076 * => djgi(xi) + Xjbi = -( digj(xj) + xfij) for z / j. (A.34) 


Dehne G : C/ —>• M n by 


GjOe) = + 



for j = 1 ,.. 


n. 


(A.35) 


Then 


diGj(x) 


digj(xj) +Xibj i^j 

0 i = j, 


(A.36) 


which together with (I A.34jl implies that BG = 0 in O. Then by Lemma IA.4I we have that 
G(x ) = a + Ax for a € M n and A = — A T € W nxn . Then 


gj(xj) = aj + ( Ax)j - (A.37) 

which we may plug into (1A.33D to deduce that 

I I 2 

u(x) = a + Ax + 7 x + (b ■ x)x — ~^~b- (A.38) 

So, every element of ker(B°) is of this form. An easy computation reveals that every function 
of this form is also in the kernel. 

□ 


Next we record a version of Korn’s inequality involving only the deviatoric part, B°. 

Lemma A. 6 . Let U C M n , n > 3, be bounded and open with Lipschitz boundary. Then there 
exists a constant G > 0 such that 

\\u\\ 2 m < C(\\u \\ 2 L 2 + ||B°u|£ 2 ) (A.39) 

for all u € H l {U). 

Proof. See Theorem 1.1 of [ 6 ]. □ 

We now prove a first result on the way to our desired Korn-type inequality. 

Lemma A.7. Let n > 3 and 

n:=n£r 1 1 (L i T)x(a,&) (A.40) 

for 0 < Li < 00 , i = 1 , ..., rt — 1 , and —00 < a < b < 00 , and define the lower boundary by 

E a = n^Tj 1 (L t T) x {a}. There exists a constant C > 0 so that 

IMIO ^ ^(IKHIo+HWa)) 


for all u G iL 1 (H). 


(A.41) 
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Proof. Suppose not. Then there exists a sequence u n G H 1 (fi) so that ||u n || 0 = 1 but 

||H5 0 u n ||o + ||' lt n||ffO(s a ) < —■ (A.42) 

Employing Lemma [A.fit we find that sup n ||u n ||^i < oo. So, up to the extraction of a subse¬ 
quence, we have that u n —*■ u weakly in H 1 and u n —>• u strongly in 77°. By (IA.42P we have 
that ID 0 u = 0 and u = 0 on E a . Then by Lemma I A. 51 we know that 

u(x) = a + Ax + 7 ® + (6 ■ x)x — b\x \ 2 /2 for y G R, a, b G M n , A = —A T G M nxn . (A.43) 

Since u vanishes on E a we must have the zeroth, first, and second order parts of the polynomial 
all vanish, and hence 

a = 0 

< Ae + ye = 0 (A.44) 

(b ■ e)e — 6/2 = 0 

for all unit vectors e G M” . The second equality in (1A.44D implies, since A is antisymmetric, 
that 

7 = 7 |e | 2 = ye • e = — Ae • e = 0 . (A.45) 

Then Ae = 0 for all unit vectors e G R n_1 , which means that A = 0 since otherwise Ae n = 1 / / 0 
implies that yi = y ■ e* = Ae n ■ e* = —e n • Ae* = 0 for i = 1 ,..., n — 1 so that y = y n e n , but then 
= Ae n • e n = 0. Choosing e = ei in the third equality in (1A.44D implies that 

6 iei = 6 2 e 2 + ... b n e n =>6 = 0 . (A.46) 

Hence u = 0 in H, which contradicts the fact that ||it || 0 = 1. 

□ 


Now we extend this to layered domains. 

Proposition A. 8 . Let n > 3, 

:= U^(LiT) x (- 6 , 0) and := n”T 1 1 (LjT) x (0, £) (A.47) 

/or 0 < Li < 00 , i = — 1 , and 0 < 6 ,l < 00 , and define the bottom boundary by 

Efc = x {— 6 } and the internal interface by E = n/j/ (LjT) x {0}. There exists a 

constant C > 0 so that 

IKII? + IMI? < c(||bV|Io + IK^-llo) ( A - 48 ) 

for all u± G with [it] = 0 along E and U- = 0 on E&. 

Proof. By Lemmas IA.61 and IA.7l we have that 

IMI? < C(\\u.\\l + ||D°u_||*) < C(||B°U_||^ + IMItfo^)) = C ll D% -llo • ( A - 49 ) 

By trace theory and the jump condition [pit]] = 0 on E, we may estimate 

II u +IIr'°(e) = ll u -lltf°(£) — C ll u -lli ^ C ||®° u -|lo • (A.50) 

Then to conclude we use Lemmas IA. 6 I and IA.7I together with this estimate: 

IKII? < c'(||u + ||2+||bV||o) < c-(||dV||S+IKII^(e)) < c/||bV||o+IK«-IIo)- ( A - 51 ) 

□ 
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A.5. Linear combinations of estimates. We now consider how to reduce a sequence of 
related estimates to a single estimate in a more useful form. 

Lemma A.9. Let k > 1 be an integer, U Cl be an open set, and consider functions Xi,Y {, Z{ : 
U —>• R for i = 0,..., k, Wj : U —> [0, oo) for i = 0,..., k, and Vq : U —> M. Suppose that 


d_ 

dt 


X 0 + Po(Yq + Wq) < 70V0 + Zq 


and for i = 1 ,..., k 


d 

dt 


i -1 


Xi + pi(Yi + Wf)< 7i ( Vq + J 2 Y 3 I + ^ 

3=0 


(A.52) 


(A.53) 


where /%, 7 * > 0 are constants for i = 0,..., k. Then there exist Aj > 0 for i = 0,..., k so that 

k k k k 


d_ 

dt 


■ Z XiXi + Z(*i + Wi) < F 0 Z Xi ^ + Z XiZi ' 

z =0 i=0 i =0 Z —0 

Moreover, A^ = l/Pk an d for i = 0,..., k — 1 

k 

a i=Pi (Pi+ 1 , ■ ■ ■, Pk, 7»+i, • • •, ik) n pj i 


(A.54) 


(A.55) 


3 =0 


where Pi is a polynomial with positive integer coefficients. 


Proof. Suppose that A* > 0 for i = 0,..., k. Multiplying (1A.53D by Ao, (I A. 531) by Aand 
summing, we find that 


d_ 

dt 


i =0 


i =0 


i =0 


k k k k k i —1 k 

Z x * x i + Z A ^ y * + Z Xi & Wi ^ £ 7 4 A 4 + 

i=1 i=0 z=l 

/c —1 Ai /c 

id) 7 aA?; + Z y * Z + Z XiZi - 

i=0 j=4+l 4=0 


i—1 

k 


(A.56) 


i=l 


We then seek to choose the A* so that 

k 


A iPi = Z Aj 7 j + 1 for i = 0 ,..., k — 1 , 
j=l+i 

XkPk = 1 - 


(A.57) 


This linear systems admits a unique solution with Aj > 0. Indeed, we may set A& = 1/Pk > 0 
and then iteratively solve for 


A,: = 




(A.58) 


for * = k — 1,..., 0. Plugging the equations (IA.57I) into (IA.56D then yields the estimate 
k k k k k 


d_ 

dt 


Z A iXi + Z + Z X ^ Wi ^ Z A * 7i + Z XiZi - 


(A.59) 


4 = 0 


4=0 4=0 


4=0 


4=0 


The estimate (IA.54I) follows from this by observing that by construction A iPi > 1 and that 
Wi > 0. The form (IA.55I) follows easily from (j A. 58 11 and elementary calculations. 

□ 
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A. 6 . Elliptic estimates. Here we consider the two-phase elliptic problem 


—pAu — (p/3 + p')SJ div u = F 2 

in H 

-§(u+)e 3 = F'l 

on E + 

- [§(«)] e 3 = -Ff 

on E_ 

N = 0 

on E_ 

u- = 0 

on Eft. 


(A.60) 


We have the following elliptic regularity result. 

Lemma A.10. Let r > 2. If F 2 € H r ~ 2 (Ll), F 3 € iL r ~ 3 ' / 2 (E), then the problem (1A.60I) admits 
a unique strong solution u € oH 1 ^) n H r (Ll). Moreover, 

" p2n 1 1111,311 (A.61) 


1 7 / II < E ll 4- UFA, 

l u llr ~ IH Hr—2 ^ lr Mr—3/2 ’ 


Proof. We refer to [33., Theorem 3.1] for the case of two-phase Stokes problem, but the proof is 
the same here. It follows by making use of the flatness of the boundaries E-|- and applying the 
standard classical one-phase elliptic theory with Dirichlet boundary condition. □ 

We let G denote a horizontal periodic slab with its boundary dG (not necessarily flat) consist¬ 
ing of two smooth pieces. We shall recall the classical regularity theory for the Stokes problem 
with Dirichlet boundary conditions on dG, 

—pAu + 'S/p = / in G 

di vu = h in G (A.62) 

u = <p on dG. 

The following records the regularity theory for this problem. 

Lemma A.11. Let r > 2. If f € H r ~ 2 (G),h €E H r ~ 1 {G),g> € H r ~ 1 / 2 (dG) be given such that 


Ik-1 ■ 

Jg Jog 


<P • A 


(A.63) 


then there exists unique u € H r (G),p € H r 1 (G)(up to constants) solving (IA.62I) . Moreover, 
IMI H r (G) + II V PIIlR.- 2 (G) < ||/||/fr-2(G) + INI/fr- 1 (G) + 11 <P 11 H r ~ !/ 2 (dG) • (A.64) 

Proof. See [T71 3T] . 

Lemma A. 12. Suppose that a > 0 and k > 0. If ( satisfies J t2 ( = 0 and 

— <tA*C + nf = ip on T 2 , 

then for r >2, 

ff IIC|| r + K||Cllr-2 ~ IMIr-2- 

Proof. We use (JA.65I) to see that 

(cr |n| 2 + n)((n) = (pin) for all n € (Lf lr L) X (Lf 1 Z). 

In particular, this means that e( 0) = 0. Since k > 0 we may estimate 


□ 

(A.65) 
(A. 66 ) 


t7 n 


C (n) < \p{n)\' 


and then argue as in Lemma 14.51 to deduce the bound 


O’ 


I r ^ 


r—2 


(A.67) 
(A. 68 ) 
(A.69) 


On the other hand, since u > 0, we may bound n 2 C(n) < |e(r )| 2 to get the estimate 

«IICII r _2<IMI r _2- (A.70) 


□ 
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